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This is a review of the most resent results from the investigation of the Ultrahigh Energy 
Cosmic Rays, particles of energy exceeding 10^* eV. After a general introduction to the 
topic and a brief review of the lower energy cosmic rays and the detection methods, the 
two most recent experiments, the High Resolution Fly's Eye (HiRes) and the Southern 
Auger Observatory are described. We then concentrate on the results from these two 
experiments on the cosmic ray energy spectrum, the chemical composition of these 
cosmic rays and on the searches for their sources. We conclude with a brief analysis of 
the controversies in these results and the projects in development and construction that 
can help solve the remaining problems with these particles. 
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I. INTRODUCTION 

Cosmic rays are defined as charged nuclei that origi- 
nate outside the solar system. Such nuclei of total en- 
ergies between one GeV to above 10^^ GeV have been 
detected. Below energies of several GeV cosmic rays 
are usually studied in terms of kinetic energy = 
Etot — mc^. In such terms the cosmic ray energy spec- 
trum extends more than 14 orders of magnitude, from 
10^ eV to above lO^o eV. 
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The exploration of cosmic rays began as a mixture 
of physics and environmental studies almost a hundred 
years ago. After the discovery of radioactivity it was no- 
ticed that between 10 and 20 ions were generated per 
cubic centimeter of air every second. The main question 
was if this ionization was a product of the natural ra- 
dioactivity of the Earth. The agent of this radioactivity 
was assumed to be 7-rays because the two other types of 
radioactive rays: a-rays (ionized He nuclei) and /3-rays 
(electrons) were easily shielded. To prove that natural ra- 
dioactivity is the culprit physicists started measurements 
of the ionization at different heights above the surface. 
Such measurements were done at the Eiffel tower. 

Just before the First World War Victor Hess started 
measuring the ionization on balloons. In 1912 he flew a 
balloon from Austria to an altitude of 5 km and to every- 
body's surprise the ionization increased by a factor of two 
rather than decrease. Werner Kohlhorster flew balloons 
to altitudes exceeding 9 km in Germany and measured 
even higher ionization level of the Hdhenstrahlung (high 
altitude radiation) as the cosmic rays were called by the 
first explorers. The term cosmic rays was put together 
by Robert Millikan, who was trying to prove that cosmic 
rays are 10 to 100 MeV 7-rays from nucleosynthesis of 
the common C and O elements. 

Kohlhorster continued his cosmic ray research during 
1930s. In collaboration with Walther Bothe he proved 
that cosmic rays can penetrate through heavy absorbers. 
Bruno Rossi shielded his detectors with one meter of lead 
and saw some cosmic rays still penetrating. Many expe- 
ditions were organized at high mountains to study the 
interactions of cosmic rays with the geomagnetic field. 
Arthur Compton organized expeditions at different ge- 
omagnetic latitudes which proved that cosmic rays are 
positively charged particles. More of them come from 
the West than from the East because the geomagnetic 
field bends positively charged particles coming from the 
West towards the surface of the Earth and those from 
the East away from it. 

Cosmic ray research was the basis for the development 
of the QED and the electromagnetic cascade theory. To- 
wards the end of the decade Pierre Auger and collabo- 
rators made several experiments at high mountain alti- 
tude where they ran in coincidence Geiger-Miiller tubes 
at large distances from each other. They concluded that 
primary cosmic rays generate showers in the atmosphere. 
Kohlhorster and Rossi ran similar experiments even ear- 
lier but of smaller dimensions. Auger estimated that the 
showers that were detected came from a primary cos- 
mic ray of energy up to 10^ GeV. The term 'shower' is 
an English translation by Patrick Blackett of the Ital- 
ian expression sciami that Rossi used in conversations 
with Beppo Occhialini. The knowledge accumulated in 
the 1930s was published in the magnificent article of 
([Rossi and Greisenl Il94l[ ) "Cosmic Ray Theory" . This 
is the beginning of the investigations of the high energy 



cosmic rays, of their energy spectrum and composition. 

Figure[T]shows the energy spectrum of cosmic rays with 
energy above 10^^ eV. Note that lower energy cosmic ray 
spectrum at Earth is affected by the magnetic fields of 
the heliosphere and the geomagnetic field. The cosmic 
ray flux as a function of energy is multiplied by E^ to em- 
phasize the spectral shape and to indicate the amount of 
energy carried by cosmic rays of different energy. This is 
a smooth power law spectrum that contains three general 
features: the cosmic ray knee above 10^^ eV, the cosmic 
ray ankle at about 3x10^® eV (3 EeV), and the cut-off 
above 3x 10^^ eV. The approximate positions of the knee 
and ankle are indicated with arrows above them. The 
cosmic ray spectrum below the knee is a power law E~" 
with spectral index a = 2.7. Above the knee the spec- 
tral index increases with Aa = 0.3. Above the ankle the 
power law spectrum becomes flatter and similar to that 
before the knee. 
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FIG. 1 Differential energy spectrum of cosmic rays of energy 
above 10^^ eV multiplied by E^. The positions of the cosmic 
rays knee and ankle are indicated with gray arrows. The ex- 
periments that contribute data to this graph are shown. The 
equivalent laboratory energy of the Large Hadron Collider is 
also shown. 



The values of the spectral indices show that below the 
knee the flux decreases by a factor of 50 when the en- 
ergy increases by an order of magnitude. Above the knee 
the decrease is by a factor of 100. Because of the de- 
crease, cosmic rays of energy above 10^^ eV are difficult 
to measure by direct experiments performed on balloons 
and satellites. The flux of such cosmic rays is about 3 
particles per hour per steradian in one square meter de- 
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tector. Particles above 10 eV can only be measured by 
air shower arrays of areas more than 10^ m^. Various 
air shower experiments obviously have different energy 
assignments that lead to the inconsistencies in the pre- 
sented spectra. 

The standard thinking in the field of cosmic rays is 
that particles of energy below and around the knee are 
accelerated at galactic astrophysical objects, mainly at 
supernova remnants and possibly at powerful binary sys- 
tems. The knee itself is probably a result of reaching the 
maximum energy of such accelerators. Particles above 
the ankle are believed to be of extragalactic origin. They 
may be accelerated at active galactic nuclei (AGN), at 
radio galaxies, in gamma-ray bursts (GRB), or in other 
powerful astrophysical systems. It is not obvious where 
the particles above the knee and below the ankle are ac- 
celerated, possibly at some special, very efficient galactic 
accelerators. 

In this article we will concentrate on the cosmic rays 
of energy above 10"'^* eV, in the lower right hand of the 
graph. The search for such high energy cosmic rays 
started in the 1950s by the MIT group led by B. Rossi. 
The first announcement of a cosmic ray shower of energy 
above 10^^ eV came from the Volcano Ranch air shower 



array in New Mexico ( (jLinslev et all 119611 ) ) that had an 



area of about 8 km^. Two years later John Linsley re- 
p orted on the detection of an event of energy 10^° eV 
r (|Linslevl . [T963h V The discoveries continued during the 
next 50 years with larger and larger arrays but the total 
world statistics is still small. 

These are the ultrahigh energy cosmic rays (UHECR) 
at least a part of which are of extragalactic origin. We 
will discuss the requirements for acceleration of such 
particles that carry more than seven orders of magni- 
tude more energy than the LHC beam and their prop- 
agation in the intergalactic space from their sources to 
us. We will introduce the UHECR detection methods 
and detectors and the results on the cosmic ray spec- 
trum and composition. We concentrate on the new re- 
sults presented by the HiRes experiment and the Auger 
Southern Observatory to which we will often refer as 
HiR es and Auger. Please consult the excellent review 



of iNagano and Watsor . 120001) for th e older experiments 
and results and that of ( Cronin . 1999f ) for the importance 
of the research in this field. So me more information could 
be found in the reviews of (iBluemer et al. . and 
( Beattv and Westerhofj . 2009h . We will conclude with a 
discussion of the remaining problems and description of 
possible future experiments. 



II. EXTENSIVE AIR SHOWERS 

Extensive Air Showers (EAS) arc the particle cascades 
following the interaction of a cosmic ray with an atom 
of the atmosphere. After this first interaction the atmo- 



sphere acts like a calorimeter of variable density with a 
vertical thickness of more than 11 interaction lengths and 
26 radiation length. 

A IQi^ eV (10 EeV, 1 EeV = 10^^ eV) proton striking 
vertically the top of the atmosphere produces at sea level 
(atmospheric thickness of 1033 g/cm^) about 3x 10^° par- 
ticles (with energy in excess of 200 KcV). 99% of these 
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FIG. 2 Heitler's schematic evolution of an electromagnetic 
cascade. At each stage of the cascade the number of particle 
is multiplied by two, either through pair creation or single 
photon bremsstrahlung. The evolution stops when individual 
particle energy fall below the critical energy, about 80 MeV 
in Air. 

are photons and electrons/positrons (referred simply as 
electrons in the following) in a ratio of about 6 to 1. 
Their energy is mostly in the range 1 to 10 MeV and they 
transport 85% of the total energy. The remaining par- 
ticles are either muons with an average energy of about 
1 GeV (carrying about 10% of the total energy), few GeV 
pious (about 4% of the total energy) and, in smaller pro- 
portions, neutrinos and baryons. The shower footprint 
(more than 1 muon per m^) on the ground extends over 
a few km^. 

The basic properties of the development of the cascade 
can be extracted from a simplified model due to Heitler. 
It des cribes the evol ution of a pure electromagnetic cas- 
cade ( (|Heitleij . [l95l l. 



A. Heitler's model of electromagnetic showers 

In his model, Heitler described the evolution of elec- 
tromagnetic cascades as a perfect binary tree (see Fig[2|). 
At each step all particles interact and produce two sec- 
ondaries of equal energy. This description assumes that 
at each step electrons split their energy in half via 
bremsstrahlung emission of a single photon while pho- 
tons produce an electron/positron pair of equal energy. 
In this simplified approach, all the processes cross sec- 
tions are taken as independent of energy and collision 
energy losses are ignored. 

The interaction step length d in the cascade is therefore 
given by the radiation length of the medium (A,. = 37 
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g/cm^ in air) as d = Xr\n2. After n steps the parti- 
cle number is Nn = 2" and their individual energy is 
Eo/Nn- This development continues until the individual 
energy drops below a critical value where the rate of en- 
ergy loss by electrons via bremsstrahlung is equal to the 
rate of energy loss by ionization. This energy is about 
E2 = 80 MeV in Air. At this point of development the 
electromagnetic cascade has reached a maximum and the 
number of particles is given by the ratio of the original 
energy to the critical one. 

Although very simplified, Heitler's model reproduces 
correctly three properties of electromagnetic cascades : 
1) The number of particles at the maximum of the cas- 
cade development is proportional to the incoming pri- 
mary cosmic ray energy : 



Nmax — Eo/E2 



(1) 



2) The evolution of the depth of maximum of the shower 
(measured in g/cm^) is logarithmic with energy: 



Xr, 



Xo + Xr HEo/EJ) 



(2) 



where Xq is the position of the start of the cascade. 
3) The rate of evolution of X^ax with energy, the elon- 
gation rate, defined as 



rflogio^^o 



(3) 



is given by the radiation length of the medium. This 
elongation rate is about 85 g/cm^ in air. 

Extensive simulations of electromagnetic cascades con- 
firm these properties although the particle number at 
maximum is overestimated by about a factor 2 to 3. 
Moreover, Heitler's model predicts a ratio of electrons to 
photons of 2 while simulations and direct cascade mea- 
surements in Air show a ratio of the order of 1/ 6th. This 
is in particular due to the facts that multiple photons are 
emitted during bremsstrahlung and that electrons lose 
energy much faster than photons do. 



B. Extension to hadronic showers 

Heitler' s model can be adapted to des cribe hadronic 
showers ( (jMatthewd . l2005t IStanevl . |2010|) ). In this case 
the relevant parameter is the hadronic interaction length 
A/ . At each step of thickness A/ In 2 it is assumed that 
hadronic interactions produce 2Nt^ charged pious and TV^r 
neutral ones. While 7r° decay immediately and feed the 
electromagnetic part of the shower, 7r+ and 7r~ interact 
further. The hadronic cascade continues to grow, feeding 
the electromagnetic part at each step, until charged pious 
reach an energy where decay is more likely than a new 
interaction. A schematic of an hadronic cascade is shown 
in Fig. [21 The interaction length and the pion multiplicity 



(3iV^) are energy independent in the model. The energy 
is equally shared by the secondary pious. For pion energy 
between 1 GeV and 10 TeV a charged multiplicity of 10 
(A'tt — 5) is an appropriate number. 
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FIG. 3 Schematic evolution of an hadronic cascade. At 
each step roughly 1 / 3rd of the energy is transferred from the 
hadronic cascade to the electromagnetic one. 

One third of the available energy goes into the electro- 
magnetic component while the remaining 2/3rd contin- 
ues as hadrons. Therefore the longer it takes for pious 
to reach the critical energy E^ (20 GeV in air, below 
which they will decay into muons) , the larger will be the 
electromagnetic component. Consequently in long devel- 
oping showers the energy of the muons from decaying 
pion will be smaller. In addition, because of the den- 
sity profile of the atmosphere, E^ is larger high above 
ground than at see level and deep showers will produce 
fewer muons. 

This positive correlation introduces a link between the 
primary cosmic ray interaction cross section with Air and 
the muon content at ground. According to those princi- 
ples primaries with higher cross sections will have a larger 
muon to electron ratio at ground. 

To obtain the number of muons in the shower one 
simply assumes that all pions decay into muons when 
they reach the critical energy. = (2iV7r)"^ where 

ric = ln(i?o In 3iV^ is the number of steps needed for 
the pions to reach E^ ■ Introducing /3 = ln2iV7r/ln3iV7r 
(0.85 for = 5) we have: 



(4) 



Unlike the electron number, the muon multiplicity does 
not grow linearly with the primary energy but at a slower 
rate. The precise value of /3 depends on the average pion 
multiplicity used. It also depends on the inelasticity of 
the hadronic interactions. Assuming that only half of the 
available energy goes into the pions at each step (rather 
than all of it as done above) would lead to /3 = 0.93. 
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Detail ed simulations Rive values o f /3 in the range 0.9 to 



0.95 ( (|Alvarez-Mufiiz et a/.l . 120021) ) 



The determination of the position of shower maximum 
is more complex in the case of hadronic shower than in 
the case of a pure electromagnetic one. The larger cross 
section and the larger multiplicity at each step will re- 
duce the value of X^ax while the energy evolution of 
those quantities will modify the rate of change of X„iax 
with energy - a quantity known as the elongation rate. In 
addition the inelasticity of the interaction will also mod- 
ify both the position of the maximum and the elongation 
rate. A proper account for the energy transfer from the 
hadronic component to the electromagnetic one at each 
step together with a correct superposition of each elec- 
tromagnetic sub-showers to compute X^ax is beyond the 
scope of a simple model but can be successfully done in 
a simulation. An approximation based on the sole evolu- 
tion of the EM cascade initiated by the first interaction 
f alls short of the fu ll simulation value by about 100 g/crn^ 
l' (|Matthewsl . l2005[) l. 

A good approximation of the elongation rate can be 
obtained when introducing the cross-section and multi- 
plicity energy dependance. Using a proton Air cross sec- 
tion of 550 mb at 10^^ eV and a rate of change of about 



50 mb per decade of energy f dUlrich et al 
obtains: 

A/ ~ 90 - 9 log (Eo/EeV) g/cm^ 



one 



(5) 



Assuming, as in (|Matthewsl . l2005l ). that the first interac- 
tion initiates 27V^ EM cascades of energy Eq/GNt^ with 
oc {Eo/PeVy^^ for the evolution of the first interac- 
tion multiplicity with energy, one can calculate the elon- 
gation rate: 



""1°= dl^ 



^10 



d{Xi \n2 + \r In [Eo/{%N^E2) 
dhgEo 



9 In 2 ~ 62 g/cm^ 



(6) 



(7) 



This result is quite robust as it only depends on 
the cross section and multiplicity evolution with en- 
e rgy. It is in good agreein ent with simulation codes 



f (|Alvarez-Mufiiz et all . I2OO2I )). 



The fast rate of the e nergy transfer in hadronic showers 
was noted long ago by ( Linslevl . ll977 ) who introduced the 
elongation rate theorem that stipulates that the elonga- 
tion rate for electromagnetic showers (-070) upper 
limit to the elongation rate of hadronic showers. This 
is of course a direct consequence of the larger hadronic 
multiplicity which increases the rate of conversion of the 
primary energy into secondary particles. 

Extension of this description to nuclear primaries can 
finally be done using the superposition model. In this 
framework the nuclear interaction of a nucleus with 



atomic number A is simply viewed as the superposition of 
the interactions of A nucleons of individual energy Eq/A. 
Showers from heavy nuclei will therefore develop higher, 
faster and with less shower to shower fiuctuations than 
showers initiated by lighter nuclei. The faster develop- 
ment implies that pious in the hadronic cascade will reach 
their critical energy (where they decay rather than inter- 
act) sooner and therefore augment the relative number 
of muons with respect to the electromagnetic component. 
From this simple assumptions one can directly see that: 
1) Shower induced by nuclei with atomic number A will 
develop higher in the atmosphere. The offset with respect 
to proton showers is simply : 



X — XP — A 

max max f 



In A 



(8) 



2) Showers initiated by nuclei with atomic number A will 
have a larger muon number : 



(9) 



3) The evolution of the primary cross section and multi- 
plicity with energy for nuclei is the same as for protons. 
Different nuclei will have identical elongation rates and 
will show up as parallel lines in an X^ax vs energy plot. 
See figure HI 

4) The fiuctuation of the position of Xmax from one 
shower to another is smaller for heavy nuclei than for 
light ones. 
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FIG. 4 Evolution of the position of Xmaa: as a function of en- 
ergy (elongation rate) for iron and proton induced air showers. 
Elongation rate of different nuclear species are with nearly 
constant slope and almost parallel to each other. Shown here 
are the results of detailed simulation performed by the Auger 
collaboration using various interaction models. 

All the above results and properties are qualitatively 
confirmed by detailed simulations. All interaction models 
share those basic principles and they all predict that iron 
showers have a smaller average Xm.ax-, less fluctuations on 
Xmax and a larger muon to electron ratio at ground than 
proton ones. In particular the offset in X^aax from iron to 
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proton showers is more than 100 g/cm^ and iron showers 
carry about 1.8 times as many muons as proton showers 
of the same energy. Of course in quantitative terms there 
are differences but all the basic trends regarding the evo- 
lution of Xjnax and with energy and atomic number 
are reproduced. This is of particular importance in the 
attempts to relate experimentally measured quantities to 
mass composition. 



C. Main features used for composition studies 

On a shower to shower basis, composition studies are 
particularly difficult because of the intrinsic shower to 
shower fluctuation of X^ax and iV^ . Those fluctuations 
come from the random nature of the interaction processes 
(in particular the position of the first interaction) and 
from the large spacing and limited sampling size of the 
detectors. Nevertheless, due to the difference in their 
cross section with Air, showers originating from different 
primaries can, at least statistically, be distinguished. 

In a real situation, where the composition evolves with 
energy, one observes changes in the elongation rate that 
are not compatible with a single species because the rate 
of change is either too large, when composition evolves 
from heavy to light (violating the elongation rate theo- 
rem) or too small going from light to heavy. 

From the superposition principle we have seen that 
distinct primaries will show up as parallel lines of con- 
stant slopes in an elongation rate plot. Detailed simu- 
lations qualitatively confirm this principle although the 
lines are neither totally parallel nor exactly of constant 
slope. Those features are model dependent as they de- 
pend on the inelasticity treatment of the cross section 
and on the leading (the fastest) particle effect together 
with the evolution of the rate of change of the cross 
section with energy which is not measured at the high- 
est energies (above 10^'' eV). Nevertheless for a detector 
that can measure in each individual cosmic ray event the 
position of Xmax with decent precision (a few tens of 
g/cm^ or less) the elongation rate plot provides informa- 
tion about the evolution of cosmic ray composition with 
energy. It is important to note that both the absolute 
value of Xmax and its rate of change with energy are not 
the same for the various interacti o n models tha t have 
been developed ((Ahn et al . 2009t Ostapchenkol 2007 : 
Pierog and WerneT 20091 )). Therefore going from the ex- 
perimental average value of Xmax at a given energy to 
an average atomic number is strongly model dependent. 

Beside the average value of Xmaj, at a given energy 
another statistical observable which distinguishes com- 
position and is less model dependent is given by the 
width of the 'K.max distribution. In a simple approach 
the shower to shower (of the same energy) fluctuations 
of y^max are dominated by the fluctuations of the first 
interaction point Xq. Xq follows an exponential distri- 



bution. As in Heitler's model, we take the 50% energy 
loss distance: A/ In 2 for the characteristic length in the 
distribution of Xq. At 10^® eV this approach gives, ac- 
cording to equation ([5]), a fluctuation of Xq for protons 
of 60 g/cm^ which is in good agreement with simula- 
tion results. For an iron nucleus of the same energy, the 
lower energy of individual nucleons and the strict ap- 
plication of the superposition principle give RMS[Xo] = 
A/(10i8eV/56)ln2/V56 ~ 14 g/cm^. This is in reason- 
able agreement with detailed simulations which give a 
value between 20 and 24 g/cm^. Leading fragment effect 
(violating the superposition principle) and fluctuation in 
subsequent interactions also play a role in the case of 
heavy nuclei. Here the fluctuations of Xq give only a 
lower bound to the fluctuations of 'K^ax ■ 

On an individual shower basis, the identification of 
the primary cosmic ray is experimentally more challeng- 
ing but not totally hopeless. Due to the fast rate of 
growth of the particle number in the cascade and to the 
large phase space available for secondaries at each in- 
teraction, the particle fluxes converge rapidly towards 
distributions that are independent of the primary par- 
ticle type. This is especially true around the maxi- 
mum development of the shower. Electromagnetic and 
muons fluxes are adequatel y described by a Gaisser-Hillas 



eqi 

C llGaisser and Hillas . Il977 )) type function (as in Eg. [T^ 



whose "age" parameter describing their stage of devel- 
opment is derived from ^max and whose normalizations 
are given by the primary energy and the muon frac- 
tion. This universali ty property was recently discussed 



m (IChou et aLl.l2005l)' and s tudied in (lApel et al 



Ciller et all 12004 iLipadl l2009t iNerling et al 
Schmidt et all l2007il2008l) 



2008; 



2006 : 



Air shower universality states that the longitudinal de- 
velopment of the electromagnetic component of nuclei- 
induced air showers can be completely described in terms 
of two parameters: the primary nucleus energy and the 
shower age. Shower universality tells us that all informa- 
tion about the primary particle can in principle be recov- 
ered from the measurements of 3 parameters. While, due 
to fluctuations, it is insufficient to measure only y^max 
and Eg , efficient separation can be achieved if the muonic 
content of the shower is also measured. Additional infor- 
mation on the first interaction cross section can also be 
retrieved by fitting an exponential to the right hand side 
(deeper side) of the Xmax distribution in fixed energy 
bins. 

Unlike the electromagnetic component of EAS muons 
reaching ground level still carry information about their 
production point along the shower axis. Because they 
mostly travel in straight lines without much scattering 
they dominate the early part of the signal at ground. 
Therefore detectors with good timing capabilities can 
construct composition sensitive parameters based on the 
signal shape either in individual detectors (rise time 
parameters) or comparing signals in several detectors 
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(asymmetries and curvature parameters) . Even when the 
absolute muon content cannot be retrieved these shape 
parameters provide valuable information characterizing 
the primary composition. 



D. Detection methods 

Above 10^^ eV the cosmic ray flux drops below a few 
tens of particle per m^ and per year. It is no longer 
possible to detect the incident particles above the atmo- 
sphere before they interact. Detectors flying in balloons 
or satellites that are less than a few m^ in size must be 
replaced by ground based instruments that cover up to 
several thousands of km^ . 

From the direct measurement of the incident particle 
properties, energy, mass, charge, etc., one must revert to 
the indirect measurement of the EAS produced by the 
interaction in Air. The atmosphere acts as a calorimeter 
and becomes part of the detection system. As this is not 
a fully controlled environment, atmospheric conditions 
must be carefully monitored and recorded along with the 
air shower data. All experiments aim at measuring, as 
accurately as possible, the primary direction (by the rel- 
ative times of the signals), the primary energy (inferred 
from the integrated signals densities), and the primary 
nature or mass (extracted from the signals shapes). 

With the exception of fluorescence light from the ni- 
trogen molecules excited along th e shower trajectory an d 
the possible microwave emission (( Gorham et al. . 2008a[ ). 
most radiations emitted from EAS are concentrated in 
the forward direction and cannot be detected far away 
from the shower axis. Hence the original, and most fre- 
quent technique, used to detect UHECR is to build an 
array of sensors (scintillators, water Cherenkov tanks, 
muon detectors, Cherenkov telescopes, ...) spread over 
a large area. When a cosmic ray event falls within the 
array boundary, the sub-sample of detectors placed near 
enough to the shower axis will observe the radiation 
reaching ground. The surface area of the array is cho- 
sen according to the incident flux, i.e. the energy range 
one wants to explore. 

Ignoring the remaining fragments of the hadronic cas- 
cade which are concentrated very near the shower core 
electrons and photons from the EM cascade, muons and 
forward beamed Cherenkov light propagate along the 
shower axis. Particles reaching ground from the EM cas- 
cade are the result of a long chain of interactions, they are 
constantly regenerated and progress in a diffusive way. 
Those observed at ground are produced in the vicinity 
(a couple of Moliere radii or a few 100 m) of the ground 
detector that measures them. Their time profiles carry 
little information on the shower development itself but 
their density gives information on the primary energy. 
This radiation is concentrated around the shower axis, 
but at the highest energies, above 10^® eV, particles can 



be observed up to a couple of kilometers away with detec- 
tors of about 10 m^ in size. Like the EM cascade, muons 
and (direct) Cherenkov light are concentrated around the 
shower axis. However, they reach ground essentially un- 
altered. Their time profile carries the memory of their 
production point along the shower axis and can be used 
to construct composition sensitive parameters. 

Fluorescence light is emitted isotropically and hence 
can be detected with appropriate telescopes tens of kilo- 
meters away from the shower axis. The light is emitted 
proportionally to the number of electrons in the EM cas- 
cade and reaches the telescopes essentially unaltered (we 
neglect the losses due to diffusion). The time profile will 
then reflect the evolution of the electromagnetic cascade 
and allows for direct measurement of composition sen- 
sitive parameters such as Xmoi- Moreover, because the 
radiation can be observed far away, a clear lateral view 
of the shower profile is possible, unlike in the case of a 
detection close to the shower axis. 



1. Air shower arrays 

Air shower arrays are networks of particles detectors. 
They cover a surface area in direct proportion to the CR 
flux in the region of the spectra one wishes to study. A 
few thousands m^ is enough for the knee region around 
10"'^^ eV, while thousands of km^ are necessary for studies 
near the spectral cutoff at energies above 10^^ eV. The 
spacing of the detector is also a function of the energy 
range of interest. For cosmic rays of energy 10^^ eV and 
above spacing is of the order of 1 km. 

The array of detectors counts the number of secondary 
particles which cross them as a function of time. They 
sample the part of the shower which reaches the ground. 
The incident cosmic ray direction and energy are mea- 
sured assuming that the shower has an axial symmetry 
in the transverse shower plane. This assumption is valid 
for zenith angles up to about 60°. At larger zenith angles 
the EM part of the cascade is largely absorbed and the 
muons start to be bent by the geomagnetic field. Above 
75°, the ground pattern shows a clear butterfly shape 
characteristic of the geomagnetic field effect. 

The pioneer work of J. Linsley at Volcano Ranch 
(' (|Linslev et al. . 196l[ )) used an array of 3 m^ scintillators 
900 meters apart covering a total surface area of about 
8 km^. It is with this detect or that the firs t event in the 
10^° eV range was detected (( Linslev . 19631) ). Scintillator 
arrays are usually made of m^ flat pieces of plastic scintil- 
lators laid on the ground and connected by cables. They 
are equally sensitive to all charged particles, thus mea- 
sure mostly the EM component of the cascade. Particu- 



larly simple to use and deploy they ha ve been quite popu 



lar for studies at the highest energies (dChiba et a/.l . ll992 : 



lEfimov et al. , 1991 ; Linslev et all . Il961 ) ) . The aperture 



of flat scintillator arrays drops quickly with zenith an- 
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gle because of the decrease of their effective surface and 
because of the absorbtion of the E.M. component. For 
accurate measurements data of scintillators array is usu- 
ally restricted to zenith angles below 45°. 

In principle the measurement of the EM cascade allows 
for a calorimetric and essentially mass independent mea- 
sure of the primary cosmic ray energy. However, detec- 
tor arrays sample the particle densities at a fixed atmo- 
spheric depth which varies from shower to shower because 
of the variations of the position of Xmax. This introduces 
a mass dependent bias in the energy estimates. In prac- 
tice the energy calibration of scintillator arrays relies on 
simulations. This has always been the major difficulty of 
the technique. 

Water Cherenkov tanks have also been successfully 
used in large cosmic ray arrays. The Haver ah Park ar- 
ray, made of Cherenkov tanks of various sizes spread 
o ver about 12 km^ t ook data for almost 20 years 
(( Lawrence et al . 199ll )). Heavy, requiring extra pure 
water with excellent protection against contamination, 
water Cherenkov datectors are not as easy to deploy as 
scintillators. However, since the Cherenkov light gener- 
ated in the water is proportional to the pathlength of the 
particle, water tanks are sensitive to both the numer- 
ous electrons and photons, and the shower muons. On 
the average, depending on the exact detector geometry, 
a muon will deposit about 10 times more light than a 
single 20 MeV electron. Because of their height, water 
tanks also offer a non zero effective surface for horizon- 
tal showers. Together with the muon sensitivities this 
extends the aperture of such arrays to nearly horizontal 
showers. 

Reconstruction of the primary particle parameters is 
based on timing for the geometry and on the distribution 
of signal densities as a function of the lateral distance to 
the shower axis for the energy. 

From the position of the different detectors and from 
the onset of the shower front signal recorded in each of 
them, one can reconstruct the shower axis and hence the 
original cosmic ray direction. Precision of one to three 
degrees are usually obtained given the large base line of 
the detector spacing (1 km). For charged cosmic rays 
this precision is sufficient given the deflection expected 
from the galactic coherent and random magnetic field. 

For the energy, the detector position are projected onto 
the plane transverse to the shower axis and a "lateral 
distributi on fun c tion" (LDF) is adjusted to the measured 
signals. (iHillad . [l970h proposed to use the signal at an 
optimal distance ropt depending on the energy range and 
the array spacing. At ropt the sum of the fluctuations 
from shower to shower and of the statistical fluctuations 
from particle counting are minimum. 

Several LDFs have been used to represent the lateral 
signal distribution. The Haverah Park experiment used 



as LDF the function: 



(10) 



for distances less than 1 km from the shower core. Here r 
is in meters, 9 is the zenith angle and /3 can be expressed 
as: 



13(6, E) ^a + bsece 



(11) 



The value of f3 is derived from Monte Carlo simulations 
and can also incorporate a slow logarithmic evolution 
with energy. Roughly a ~ 3.5 and b ~ —1.2 for a verti- 
cal value of f3 under 3. At larger distances (and higher 
energies), this function has to be modified to take into ac- 
count a change in the rate at which the densities decrease 
with distance. This is due to the increasing dominance of 
muons over electrons at large core distance, since muons 
have a flatter LDF. A more complicated for m was used 



by the AGASA group f llYoshida et aLl . ll995l )). However 



the principle remains the same. On figure [5] an exam- 
ple, taken from the Auger collaboration public event dis- 
play Q of the footprint of an air shower on the ground 
together with the reconstruciton of the LDF is shown. 




J 5 1 500 2000 
Distance to axis [m 



FIG. 5 Example of detection using a surface array. The 
upper right insert shows the whole auger surface array and 
the footprint of the shower, each dot represented a detec- 
tor the spacing between them is 1.5 km. The lower insert 
shows details of this footprint with the estimated contours 
of the particle density levels. The curve represent the ad- 
justed LDF (lateral distribution function) and the red point 
the measured densities as a function of the distance to the 
shower core, (extracted from the public event display of the 
Auger collaboration). 

Once the attenuation of the signal due to the zenith 
angle is accounted for, an estimator of the energy is ob- 
tained from the corrected density at ropt in the form : 

E = kSiropt,erefr 



(12) 



http:/ /auger. colostate.edu/ED/ 
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where a is a parameter of order 1. 



2. Cherenkov light 



To reconstruct the primary parameters, a minimum 
of three detectors with signal is necessary. The spacing 
between those detectors will determine the array energy 
threshold. For a vertical shower the 500 m spacing of 
the trigger stations in Haverah Park corresponds to a 
threshold of a few 10^^ eV, while the 1.5 km separation 
of the Auger Observatory stations gives a few 10^* eV. 

Ground arrays do not have a direct access to the posi- 
tion of the shower maximum and this is a strong limita- 
tion to this technique for primary identification. Muon 
counting can be done with buried detectors or, in favor- 
able conditions, when the EM to muon signal ratio is not 
too large, by counting muon spikes in the recorded traces 
of water Cherenkov tanks. Additionally, when again the 
EM to muon signal ratio is not too large, the early part 
of the traces is dominated by the muon signal and its 
time evolution carries information on the position of the 
shower maximum. This is always less sensitive than a 
direct measurement of X^ax as done by the fluorescence 
technique. Only a combination of both measurements, 
Xmax and Nfj_ can give a shower by shower composition 
indication. 

Alternative techniques trying to exploit the emission 
of EAS in the radio band have also been explored. Be- 
tween 1967 and 1973 extensive studies took place in the 
10-100 MHz band. However the technique was judged 
unworthy in particular due to the strong beaming of the 
emission in the forward direction an d to the poo r signal 
to noise ratio achieved at the time (|Allanl . ll977h . With 
the progress in fast digital electronics and low noise am- 
plifiers the i nterest for the t echnique was revived. The se 
new efforts f (|Ardouin et ai . , I2OO6I: iHaungs et all l2009l) ). 
aim at replacing ground array detectors with radio an- 
tennas which are both less expensive and easy to deploy. 
In addition, the radio signal which propagates essentially 
non altered from its source to the detector carries infor- 
mation on the shower evolution. Important progress has 
been made and the radio signal in the VHP band has 
been showed to be dominated by the geo-synchrotron 
emission of the electrons and positrons of the EM cas- 
cade. However, detection of transient signal in those fre- 
quencies is still a challenging task, even more so given 
the very tight lateral distribution of the radio signal. 
Recent measurements confirm the strong concentration 
of the signal in the forward region with an exponential 
decrease from the core with a characteristic distance of 
order 150-200 m. However they demonstrated the pos- 
sibility to reconstruct the CR direction with reasonable 
accuracy (( Revenu et al . 2009l )l. Progress regarding this 



Ac cording to Brennan and Chudakov f ()Brennan et al. 



technique and its exploitation at the highest energies are 
expected in the co ming years as import ant R&D efforts 
are being pursued (I van den Berg . 20091) '). 



19581 : IChudakov et all the Cherenkov light emis- 

sion from the charged particle component of an air shower 
can provide an integrated measurement of the longitu- 
dinal development. The Cherenkov intensity is propor- 
tional to the primary energy, while the slope of the lateral 
distribution is related to the depth of maximum shower 
development. Thus, if one samples the Cherenkov lat- 
eral distribution, i.e. the photon density as a function 
of the distance from the air shower core, it is possible to 
estimate both the primary energy and composition. 

Prom air shower simu lations, it was shown 
(( Patterson and Hillas . 19831 )) that at distances larger 
than about ^ 150 m from the shower core, the density 
of Cherenkov light is proportional to the primary energy 
but essentially independent of its nature. The light 
profile close to the core is sensitive to the depth of 
penetration of the shower in the atmosphere which 
correlates with the primary cross section. 

Experimental setups exploiting the Cherenkov light 
for EAS are usu ally associated to standar d particle 
detec t or arrays ((Araueros et al. [2OO O: Ca ssidv et al 



1997 



1999; 



.Chernov et all 120061: I Disckinson et al. . 

Efimov et all Il99ll: ISwordv et all \200lh ). Cherenkov 
light is also used in CR observation at other energies. 
For a complete overvie w and history of C herenkov detec- 



tion of cosmic rays se e (iLidvanskyI, l2005h 



In dCassidv et all Il997l) . (|Chernov et all I2OO6I ) 



and ( Araueros et aLl . l2000l ) the experimental setups con- 
sist of open photomultipliers fitted with Winston cones 
and looking upward in the sky. The largest array com- 
posed of 150 of such sensors distributed about every 
40 m was installed on the Ply's Eye site in Dugway 
(Utah, USA) to gether with the CASA and CASA-MIA 
(^ (jSorione et all . Il994) ) detector arrays. Near the core, 
the lateral distributio n of Ch erenkov light was shown to 
be exponential as in (IPatterson and H illas. 1983). The 
C ASA-BLANC A group f dPowler et oi.l . l200lh ) used a two 
component function which matches both their real and 
simulated data. The function is exponential in the range 
30m-120m from the shower core and a power law from 
120m-350m. It has three parameters: a normalization 
C120, the exponential "inner slope" s, and the power law 
index /3: 



r( '^120 e<i20m-.)^ 30m < r < 120m 

^' \ C120 (r/120m)-'3, 120m< r < 350m ^ ' 

The primary energy depends primarily on C120, the 
Cherenkov intensity 120 m from the core. Detailed simu- 
lation of the shower and of the detector can be used to de- 
rive the relation between those two quantities. Hadronic 
models predict that C120 grows approximately as E^'^'^ , 
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because in hadronic cascade the fraction of primary en- 
ergy directed into the electromagnetic component in- 
creases with energy. 

Similarly the slope of the exponential can be related 
to the shower maximum using simulatio ns. The relation 



between the two is esse ntially linear ( ([Argueros et al. 
2nnnl : iFowler et oi.l . l200lh V 



The low duty cycle (Cherenkov detector can only be 
operated on clear dark nights), the short core distance 
up to which the inner slope parameter can be used to es- 
timate Xmax and consequently the small spacing within 
units made this technique inappropriate to study EAS 
beyond an energy of about 10^^ eV. The success of the 
fluorescence detection technique contributed to the de- 
cline in the interest for this technique at the highest en- 
ergies. 



3. Fluorescent light 

The charged secondary particles in EAS produce ul- 
traviolet light through nitrogen fluorescence. Nitrogen 
molecules, excited by a passing shower, emit photons 
isotropically into several spectral bands between 300 and 
420 nm. As discussed above, a much larger fraction of 
UV light is emitted as Cherenkov photons. But this emis- 
sion is strongly beamed along the shower axis and usually 
considered as a background to fluorescence detection. 




FIG. 6 Sketch of the detection principles of a fluorescence 
detector. The fluorescence light emitted by the air shower 
is coUected on a large mirror and focussed onto a camera 
composed of photo-multipliers (Auger collaboration). 

The first fluorescence detector assembled for UHECR 
detection was laid down by Greisen and his te am in the 
mid 60's ( (iBunnej . Il967t bunner et all Il967t )). SmaU 



mirrors and the atmospheric conditions did not allow to 
record signals from EAS. Detectors were built in the late 
70 's by a group of the University of Utah and tested 



at th e Volcano Ranch ground array ( (jBergeson et al. 



fl Hara et al . 1970f) ). Later on, a fully functional detec- 
tor was ins talled at Dugway (Utah ) under the name of 
Fly's Eye (( Baltrusaitis et adilQSSi) ). It took data from 
1981 until 1993 and fully demonstrated the extraordi- 
nary potential of the technique . The highest energy 
shower ever detected (320 EeV) was observed by this de- 
tector. An updated version of this ins trument, the High- 
Resolution Fly's Eye, or HiRes ( (Bover et a/.l . l2002f )). ran 
on this same site from 1997 until 2006. 

The fluorescence yield is 4 photons per electron per me- 
ter at ground level pressure. Under clear moonless night 
conditions, using square-meter scale telescopes and sen- 
sitive photodetectors, the UV emission from the highest 
energy air showers can be observed at distances in excess 
of 20 km from the shower axis. This represents about two 
attenuation lengths in a standard desert atmosphere at 
ground level. Such a large aperture, instrumented from 
a single site, made this technique a very attractive al- 
ternative to ground arrays despite a duty cycle of about 
10%. 

Fluorescence photons reach the telescopes in a direct 
line from their source. Thus the collected image reflects 
exactly the development of the EM cascade (see flgure[6]). 
From the fluorescence proflle it is in principle straight- 
forward to obtain the position of the shower maximum 
and a calorimetric estimate of the primary energy. In 
practice a number of corrections must be made to ac- 
count for the scattering and the absorption of the fluo- 
rescence light. Also pollution from other sources such as 
the Cherenkov component which can be emitted directly, 
or diffused by the atmosphere into the telescope, must be 
carefully evaluated and accounted for. A constant moni- 
toring of the atmosphere and of its optical quality is nec- 
essary together with a precise knowledge of the shower 
geometry for a careful account for those corrections. 




19771) ) while the flrst detection of fluorescence light from 
UHECR was made by Tanahashi and his collaborators 



FIG. 7 Geometry of the detection of an air sho wer by a fluo- 
rescence telescope, from jKuempel et aZ.I . [20081 ) . 

The shower geometry as viewed from a fluorescence 
telescope is depicted in Fig. [71 It is defined by the shower 
detector plane (SDP), the distance of closest approach 
i?p, the time io along the shower axis at the distance 
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of closest approach and the angle xo within the SDP 
between the ground plane and the shower axis. This 
geometry is usually reconstructed in two steps. First, the 
shower-detector-plane (SDP) is determined. Next, the 
arrival time U of the signal in each pixel in the direction 
of SDP Yi is used to deterni i ne yn , Rp and to from the 



Gaiss er-Hillas form and a three -parameter Gaussian in 
age (' (|Abu-Zavvad et al. . 200 1[ )). Alternatively, ana- 



equation f ()Baltrusaitis et ami985D ) 



«), -ftp 



Rp 
c 



ti —to + — tan f 



(14) 



One important property of this equation is that un- 
less the angular velocity in the camera and its rate of 
change can be measured there is a degeneracy between 
the impact parameter Rp and the angle xo- This de- 
generacy leads to poor pointing resolution - the three 
parameters defining t he shower g eome try cannot be con- 
strained accurately ( dSommers . 1995f) V The situation 
can be improved using fast electronics to achieve a good 
precision on ti and for those showers with sufficient track 
length in the camera (over about 10°). This was first 
used by the HiRes cohabora tion with the HiRes-II de- 
tector (( Abbasi et aLl . lioOSah ). Alternatively the HiRes 



collaboration also developed a profile constrained time fit 
(PCF) for t he part of its detecto r not equipped with fast 
electronic f( Abbasi et al. . 2005al )). Nevertheless, in both 
cases the geometrical resolution remains at a few degrees. 
The best option to resolve this ambiguity is to improve 
the measurements. This can be done in two ways: 

A) Using a second telescope viewing the shower from a 
different position. The intersection of the two SPD will 
constrain the geometry of the shower axis to within a 
fraction of a degree. This is called a stereo reconstruc- 
tion and is the technique used by the HiRes detector. 

B) Constrain the to parameter by a direct measurement 
of the time of arrival of the shower at the ground. This is 
the hybrid technique used by the Auger detector. Again 
the geometry can then be constrained to within a fraction 
of a degree. 

Once the geometry has been determined, the fiuores- 
cence technique is the most appropriate way to measure 
the energy of the incident cosmic ray. The amount of 
fluorescence light emitted along the shower axis is pro- 
portional to the number of electrons in the shower. The 
EAS has a longitudinal development usually parameter- 
ized by the 4 parameter Gaisser-Hillas function giving 
the size Np oi the shower as a functio n of the atmospheric 
depth X C dGaisser and HillaslllQTi '): 



X-Xq 

X-ma-x — X(] 



-Xo)/\ 



=.(X„„„-X)/A 



(15) 

The total energy of the shower is proportional to the 
integral of this function, knowing that the average energy 
loss per particle is 2.2 MeV/g cm^^. 

The Pierre Auger fluorescence reconstruction uses 
this formula while the HiRes group has used both the 



lyti c shower theor y led to yet another form popularized 
bv (iGreisenl . llQSdl ). In a recent study f ([Matthews et al . 
\20iw ) it was shown that the introduction of the pro- 
file full width half-maximum { fwhm) and its asymmetry 
(defined by the ratio of the left-width at maximum to 
the fwhm) could unify the parameterization of all three 
profile functions. Greisen and Gaisser-Hillas profiles are 
shown to be essentially identical while gaussian in age 
profile only differ at the very early and very late devel- 
opment stages of the cascade. 

Beside the corrections arising from the experimental 
conditions discussed above, the energy transported by 
the neutral particles (neutrinos), the hadrons interacting 
with nuclei (whose energy is not converted into fluores- 
cence) and penetrating muons, whose energy is mostly 
dumped into the Earth, must also be accounted for to 
estimate properly the primary CR energy. This missing 
energy correction is calculated using detailed simulations 
and varies with energy, composition and the interaction 
model used. It is about 20% at lO^^ eV for iron (10% 
for proton) and about 12% at 10^*^ eV (6% for protons). 
Variations from one model to another are of about 50% 
( (jPierog and Werneil . l2007h ). 

Despite the fact that fluorescence measurements give 
direct experimental access to the position of X^ax the 
separation of hadronic primaries according to their mass 
cannot be done on a shower by shower basis because of 
the intrinsic fluctuation of this parameter. One must look 
for statistical means of studying the chemical composi- 
tion and / or use additional information such as the muon 
content that can be provided by particle detectors as in 
the hybrid detection system. 

Statistical methods relying on the measured fluores- 
cence profile only are based on the elongation rate plot, 
or the RMS{Xmax) plot. In the former, one calculates the 
average value of X^ax from a set of showers of the same 
energy and plot it as a function of energy. In the latter, 
it is the width of the distribution of Xmax of shower of 
the same energy that is plotted against energy. Those 
measurements are discussed in section IVIll 



III. GALACTIC COSMIC RAYS 

A. Origin of the galactic cosmic rays 

Galactic cosmic rays are believed to be acceler- 
ate d at supernova remnants. Thi s idea was justified 
by (jGinzburg and Svrovatskiil |l96J) through simple and 
powerful arguments based on the energetics of supernova 
remnants. If only 5% to 10% of the kinetic energy of su- 
pernova remnants is converted to accelerated cosmic rays 
this would provide the energy of all galactic cosmic rays. 

Supernova remnants are attractive candidates for cos- 
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mic ray acceleration because they have higher magnetic 
fields than the average interstellar medium. They are 
also large and live long enough to carry the acceleration 
process to high energy. The acceleration mechanism is 
believed to be stochastic acceleration at supernova blast 
shocks. 

The idea of stochastic particle acc eleration was first de- 
veloped by E. Fermi who proposed (^ (lFermil . [l949[) l to use 
the charged particle interactions with interstellar clouds 
to accelerate cosmic rays. 

The shock ahead of the expanding supernova remnant 
is formed because the expansion velocity of the rem- 
nant is much higher than the sound velocity of the in- 
terstellar medium. Shock acceleration is much faster 
than the original Fermi acceleration mechanism. The 
energy gain is proportional to /? (first-order accelera- 
tion) rather than to (second-order (Fermi) acceler- 
ation) where /3 is the velocity of the magnetic cloud or 
the blast shock velocity in terms of c. In addition, the su- 
pernova shock velocity is much higher than the average 
velocity of molecular clouds. As a result shock accel- 
eration is orders of magnitude more efficient, and cor- 
respondingly much faster. The shock acceleration sce- 



nario was suggested in the late 1970s ( (jAxford et al. 



19771 : iBeiil . Il978t iKrvmskv*. 119771)1 and is under contin- 
uous developinent ((.Blandford and Eichleil.ll987 : Drurv . 
19831: IJokipii Il987t Ijones and Ellisonl Il99l[ )). The pre- 
diction is for a flat, E~^ cosmic ray spectrum in ac- 
celeration in non relativistic shocks and for a steeper 
g2.2-2.3 spectrum at acceleration in highly relativistiv 



shocks f (|Achterberg et a/.l . liooTbh ). 



The maximum energy that a charged particle could 
achieve is then expressed as a function of the shock veloc- 
ity and extension and the value of the average magnetic 
field as 



Eraax = f^ZeBrs, 



(16) 



where /? is the shock velocity in terms of the speed 
of light, rs is the shock radius, and Ze is the parti- 
cle charge. Equation (ITB)) is valid during the period of 
the free expansion of the supernova remnant when the 
shock velocity is constant. During the Taylor-Sedov 
phase, when the shock has collected enough interstel- 
lar matter to start slowing down, the maximum energy 
starts decreasing as the radius is only proportional to the 
time t o the power of 2/ 5. Detailed more recent calcula- 
tions ( ( Berezhkol 19961 )) derive maximum energy values 
close to 5 X 10^ Ge V and even higher in some cases 



f (|Ptuskin et al 



2010[ )). An important component of the 
expression for Emax is its dependence on the particle 
charge Z. It means that a fully ionized heavy nucleus 
of charge Z could achieve Z times higher energy than a 
proton. 

Since cosmic rays scatter in the galactic magnetic fields 
we cannot observe them coming from particular sources. 
The only way we can study their acceleration sites is by 



observing the neutral particles, gamma-rays and neutri- 
nos, generated by their interactions during acceleration. 
There are two epochs in supernova remnant evolution 
when one can expect 7-ray and neutrino emission. One 
of them is shortly after the supernova explosion, when 
the density of the expanding supernova envelope is very 
high and thus contains enough of a target for hadronic 
interactions. The emission will continue for about 2 to 10 
years (depending on the mass distribution and expansion 
velocity of the SNR) until proton energy loss on inelastic 
interactions becomes dominated by the adiabatic loss due 
to the SNR expansion. The 7-ray emission will fade for a 
long time, until the SNR reaches the Sedov phase, when 
most of the galactic cosmic rays are accelerated. Since 
this phase lasts for more than 1,000 years there should be 
many supernova remnants that are gamma-ray sources. 
The modern expectations of the 7-ray e mission of ma- 



ture supernova remnants was developed by (jDrurv et al 



19941 ). The assumption is that cosmic rays at the source 
have a much fiatter spectrum than the one observed at 
Earth as acceleration models suggest. As an example 



of th e expectations from a concrete SNR (jDrurv et al 



1994[ ) apply the calculation to the Tycho (1572) super- 
nova remnant which should be close to the Sedov phase. 
One can take the average supernova energy and density 
from different estimates i?sN = 4.5 ±2.5 x 10^" ergs, mat- 
ter density rti = 0.7 ± 0.4 and estimate the 7-ray flux for 
conversion efficiency 9 = 0.2 and distance d — 2.25 ±0.25 
kpc. The expected flux is 

-12 f Ej 



F{> E^) = ~ 1.2x10" 



VTeV 



— 9 —1 

cm ^ 



(17) 



The detection of such a flux is easily within the capabili- 
ties of the last generation of 7-ray Cherenkov telescopes. 

Figure [8l compares th e positions of supernova remnants 
from the (Green, 20091) catalog with the positions of GeV 
g amma-ray sources from the Fermi/LAT observations 
f (|Abdo etali l2009| )) and of TeV sourc es from the TeV- 
Cat catalog http:/ /tevcat.u chicago. edu\ . One can clearly 
see several coincidences. There are others that are more 
difficult to find by eye because there are so many su- 
pernova remnants close to the galactic center at very 
low galactic latitude. Many of the TeV sources come 
from the HESS survey of the galactic plane of galac- 
ti c longitudes from -30° to 30° and latitudes below 3° 



f( Aharonian et al . 2006bf )). The names of some of the 



SNR that emit TeV 7-rays are indicated in the figure 
whenever possible. There are though no gamma-rays 
coming from the Tycho supernova remnant. 

The number of direct coincidences of the supernova 
remnants locations with the directions of the gamma- 
ray sources is relatively small. What is the fraction of 
the supernova remnants that are 7-ray sources and thus 
are cosmic ray accelerators? The small fraction of 7- 
ray producing SNRs creates doubts that galactic cos- 
mic rays are generated at these objects. This may be 
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FIG. 8 Comparison of the positions of supernova remnants (x's) with GeV (Fermi/LAT: triangles) and TeV (circles) gamma-ray 
sources. 



true, but the HESS group put together an alternative 
explan ation of this effec t in th eir study of the galactic 
ridge (( Aharonian et al. . 2006al )l. Hadronic gamma-ray 
production is only possible when the matter density of 
the medium is much higher than 1 cm~'^. A very likely 
gamma-ray production site is the location of dense clouds 
of matter close to an acceleration site of cosmic rays. 
HESS observed that the peaks of the 7-ray emission from 
the region of the galactic center ridge, after subtraction 
of known sources, coincide with the positions of molec- 
ular clouds with a matter density of hundreds per cm'^. 
The total amount of mass in these clouds is 2-4 x 10^ 
solar masses. In addition, the energy spectrum of the 
7-rays is about E"^ "^ which is likely to happen close to 
the cosmic rays acceleration site. This observation may 
explain the fact that many sources of TeV 7-rays do not 
exactly coincide with the positions of SNR where the cos- 
mic rays that produce them are accelerated, rather with 
close by molecular clouds. Higher energy cosmic rays 
diffuse faster away from their sources. For this reason 
it is possible that a molecular cloud could be a source 
of TeV 7-rays before it becomes a strong source of GeV 
gamma-rays. 

TeV 7-rays have been detected from the Crab nebula 
and the supernova remnants SN1006, Cas A, RX J1713.7- 
39466, RX J0852.0-4622, W28, W48, RCW 86 and oth- 
ers. The Crab nebula is the standard candle in TeV 7-ray 
astronomy; it has a steady flux which is used to measure 
the fluxes of other sources. The models that explain best 
its gamma-ray emission do not involve hadronic inter- 
actions. They are electromagnetic models that rely on 
electron acceleration and the inverse Compton process. 



B. Energy spectrum and composition at the knee 

The energy range in which the cosmic ray spectrum 
changes its slope is called 'the knee'. Its existence was 
fi rst suggested by the Moscow S tate University group 
( (|Kulikov and Khristiansenl . 119581 )) on the basis of their 
air shower data. Many groups have studied the knee re- 



gion and the change of the cosmic ray spectrum is well es- 
tablished. Up to an energy of 10^ GeV the spectrum of all 
cosmic ray nuclei is a power law with differential spectral 
index a of 2.70-2.75. The spectral index increases by Aa 
of about 0.3 above the kn ee. A flattening of the spectra l 
index has been detected (( Ahn et al. . 20l(k Panov et al. 
201l[ )) just before the knee. 



There is no lack of theo retical ideas about the origin 
of the knee. (jPetersl . [T959h suggested that the knee is a 
rigidity-dependent effect. Rigidity R is the ratio of the 
particle momentum to its charge. It could be related 
to the maximum rigidity that can be achieved in accel- 
eration processes or to rigidity-dependent escape of the 
cosmic rays from the Galaxy. Rigidity-dependent effect 
is an attractive idea. We know that heavy charged nuclei 
can achieve Z times higher energy at acceleration. So 
a natural assumption could be that at the approach to 
the knee cosmic ray sources can not accelerate protons 
to higher energy. Then the next nucleus. He, takes over 
and the process continues in order of charge until at some 
higher energy galactic cosmic rays contain only iron nu- 
clei. Mostly the common nuclei of H, He, C, O, Si, Mg 
and Fe are represented in the cosmic ray spectrum. 

Figure [9] shows a very simple flux model with rigidity- 
dependence. It uses the spectrum and composition at 
low energy and extends it to high energy with exponen- 
tial cutoff in rigidity at 10^ GV. The thin lines show the 
contribution of different nuclear groups to the all-particle 
spectrum. The proton spectrum turns over at 10^ GeV 
and those of heavier nuclei turn over at energies of Z x 10^ 
GeV. At energies above 10® GeV there are only heavy, 
high Z, nuclei in the cosmic ray flux. The end of the 
modeled spectrum is where the Fe component is also ex- 
ponentially cut off. 

This simple model agr ees pretty wel l with the measure- 
ment s of the Kascade (( Antoni et al . 20051 )) and Tibet 
III ( (jAmenomori et all \20oS\ ) air shower arrays. The 
normalization of these experiments is slightly different 
which affects both the magnitude of the flux when it is 
multiplied by E^-"^^ and the position of the knee. The 
analysis of air shower data depends on the hadronic inter- 
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FIG. 9 A simple model of the knee which extends the 
low energy spectrum and composition to high energy with 
an exponential cutoff at 10^ GeV. The model is normal- 
ized to the all particle fl ux m easured by the Kascade ex- 
periment f ijAntoni et aZ.I . l2005l 'l) which is shown with full 
s quares. The measurem ents of the Tibet III experiment 
( (lAmenomori et all\200§D ] is shown with empty squares. 



action models used in the simulations. The dependence 
is stronger for the Kascade experiment which is located 
much lower in the atmosphere. Tibet III is close to the 
depth of shower maximum X^ax where the ratio between 
shower electrons and muons is at its maximum. This ra- 
tio decreases with the atmospheric depth but predictions 
depend strongly on the hadronic interaction model. The 
normalization of both spectra, however, depends on the 
cosmic ray composition. 

The cosmic ray composition estimated from air shower 
data is usually presented as the average value of the loga- 
rithm of the primary particle mass (In A) . Different com- 
position estimates are not in very good agreement. As an 
illustration Fig. [TU] presents th e results from the a naly- 
ses of data fr om the Kascade ((lAntoni et al . 2005 )) and 
EAS-TOP fl'Agli etta etail . |2004|) ) experiments. Both 
composition results come from the ratio of the shower 
muon density at predefined distances from the shower 
core as a function of the total number of electrons in the 
shower, N^.. These two measurements are in a fairly good 
agreement. 

Figure [10] shows that the composition becomes signif- 
icantly heavier with increasing energy. It is fully consis- 
tent with the rigidity dependent idea. The simple com- 
position model, however, does not describe the data well. 
It predicts heavier composition at 10^ GeV and lighter 
composition at 10^ GeV. A better model would require a 
different low energy composition and possibly lower max- 
imum rigidity. 

Although the majority of the experiments measure a 
cosmic ray composition that becomes heavier between 
5x10® and 10^ GeV, it is difficult to draw a definite con- 
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FIG. 10 Results from studies of the cosmic ray composition 
in the region of the knee compared to the predictions of the 
simple model presented in Fig. [S] 



elusion about the exact changes of the cosmic ray spec- 
trum and composition at the knee. All experiments agree 
that the cosmic ray spectrum steepens above 10® GeV. 
The exact position of the spectral change and the width 
of the transition region are not yet well determined. The 
composition studies, both with surface air shower arrays 
and with optical detectors, indicate a change in the av- 
erage mass of the cosmic ray nuclei after the steepening 
of the spectrum, once again with large uncertainty in 
the energy range and shape. All these numerous data 
sets are consistent with rigidity dependent effects, either 
in the cosmic ray acceleration or in their propagation. 
This second scenario assumes that lower rigidity nuclei 
are contained in the Galaxy longer. 



IV. ORIGIN OF COSMIC RAYS UP TO 10^'' EV 

The question of how to accelerate cosmic rays up to 
10^" eV has been pending since their very first observa- 
tion in the 1960's. More than thirty years later, in the 
mid 90's, the data collected by the AGASA and HiRes ex- 
periments generated a profusion of ideas. Some of them 
aimed at an explaination of the possible absence of the 
GZK cutoff and the lack of visible astrophysical sources. 
All ideas tried to find a solution to the basic problem of 
how to transfer efficiently a macroscopic amount of en- 
ergy, of the order of 20 Joules, to a microscopic particle. 

To circumvent this difficulty one of the main axis of re- 
search was the mere suppression of the accelerator itself. 
Particles are not accelerated as such but directly pro- 
duced, via the decay of some supermassive relic of the 
Big Bang, or by the collapse of topological defects, with 
energies in excess of 100 EeV. While attractive from the- 
oretical point of view, these models had the disadvantage 
of replacing the acceleration problem with the question 
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of the nature and existence of such top-down sources. 

With the observational facts cohected by HiRes and 
Auger in the past decade the situation has been greatly 
clarified. A cutoff in the high energy end of the spectrum 
is clearly visible and the limits on the fraction of photons 
and on the flux of high energy neutrinos have strongly re- 
duced the interest in the top-down models. On the other 
hand the possibility of a dominant iron component in the 
very end of the energy spectrum decreases by a factor 26 
the hard conditions placed on " standard" bottom- up cos- 
mic accelerators to reach the 100 EeV barrier. 

Nevertheless, after many decades of investigation the 
problem has not been solved, and even the extragalactic 
natu re of the sourc e s above 3 EeV has been challenged 
ffWi ck et adl2004l) . ()Calvez et all MlO)). In the follow- 
ing we briefly reviews some of the necessary conditions 
for the acceleration of UHECR at astrophysical sites and 
enumerate some possible candidates. We also briefly re- 
view the main characteristics of the top-down models. 
More deta ils on this subject can be found in the recent 
review of (jKotera and Olintol boill) , now m press. 



A. Possible acceleration sites 

Acceleration at astrophysical sites may occur princi- 
pally through two distinct mechanisms: diffusive shock 
acceleration, based on the Fermi mechanism and one 
shot acceleration in very high electric fleld generated 
by rapidly rotating compact magnetized objects such as 
young neutron stars. 

Diffusive acceleration takes place near shock waves and 
rely on the repeated scattering of charged particles on 
magnetic irregularities back and forth across the shock. 
In the case of non relativistic shock velocities the en- 
ergy gain at each crossing is of the order of AE ~ E. 
To reach energies above 1 EeV large acceleration regions 
and/or highly relativistic blast waves are necessary. In 
the case of relativistic shock the energy gain reaches 
FjE where Ts is the shock bulk Lorentz factor. Such 
g ain appears, however, to be limited to the first crossing 
( (jAchterberg et a/.l . lioOlal ) ) . 

One of the principal advantages of the diffusive shock 
acceleration mechanism is that it naturally provides a 
power low spectrum whose predicted index 7 is within 
the range of the experimental measurements. Depending 
on the exact geometry of the shock and on its relativistic 
nature, the combination of the energy gain per crossing 
and of the escape probability leads to a power law index 
of exactly 2 for the case of a strong non relativistic shock 
in an ideal gas and to indexes between 2.1 and 2.4 for 
relativistic shocks. 

( Hillasl . lT98i ) summarized the conditions on potential 
acceleration sites using a relation between the maximum 
energy of a particle of charge Ze and the size and strength 
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FIG. 11 Hillas plot for candidate acceleration sites, relating 
their size and magnetic field strength. To accelerate a given 
particle species above 100 EeV objects must lie above the 
corresponding lines. 



of the magnetic field of the site: 
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where /3 represents the velocity of the accelerating shock 
wave or the efficiency of the acceleratoi0. We show in 
Fig.[Tl]the now famous "Hillas plot" illustrating this con- 
dition. 

Looking at the Hillas diagram one sees that only a few 
astrophysical sources satisfy this necessary, but not suffi- 
cient, condition. Among the possible candidates are neu- 
tron stars and other similar compact objects, large-scale 
shocks due to merging galaxies or clusters of galaxies, 
the core and jets of Active Galactic Nuclei (AGN), hot 
spots of Fanaroff- Riley class II (FR-II) radio galaxies and 
processes associated with Gamma Ray Bursts (GRB). 

AGN have long been considered as potential sites 
where energetic particle p r oduct i on might ta ke place 
r (|Ginzburg and Svrovatskl. Il964l: llliiia^ . 1 19841) 1. AGN 
jets have dimensions of the order of a fraction of a par- 
s ec with magnet i c field of the order of a few Gauss 
( (|Halzen and Zasl . Il997l )). These parameters could in 



^ In the case of a relativistic shock the bulk Lorentz factor 
enters the right hand side of this equation ( jAchterberg et al\ , 
l2001a|)) 



16 



principle lead to a maximum energy for protons of a 
few tens of EeV. Similarly, AGN cores with a magnetic 
field of order 10^ G and size of a few 10 ~^ parsec can 
reach about the same energy. However those maxima, al- 
ready marginally consistent with acceleration up to 100 
EeV, are unlikely to be achieved under realistic condi- 
tions. The high radiation field around the central en- 
gine of an AGN is likely to interact with the accelerated 
protons while energy losses due to synchrotron radia- 
tion, Compton processes, and adiabatic losses will also 
take place. The situation is worse for nuclei that will 
photodisintegrate even faster. Such processes may lead 
t o a maximum energy o f only a small fraction of EeV 
( (jBhattacharjee and Sig 1, [2000)). To get around this 
problem, the acceleration site must be away from the ac- 
tive center and in a region with a lower radiation density 
as in the terminal shock sites of the jets, a requirement 
possibly fulfilled by FR-II galaxies. 

The link with GRB and UHECR acceleration was 



i nitially made by (IWaxmanL ^99^, (jVietril . Il995h . and 



( Milgrom and UsovT 1995h . who pointed out that the ob- 
served at the time cosmic ray flux beyond 100 EeV (now 
estimated to be lower by a factor 3 to 10 after the mea- 
surements of Auger and HiRes) is consistent with a sce- 
nario in which these particle are produced in GRB's pro- 
vided that each burst produces similar energies in gamma 
rays and in high energy cosmic rays. From a phenomeno- 
logical point of view, based on the gamma-ray observa- 
tions, bursts can be described by the product of the dis- 
sipation of the kinetic energy of a relativistic expanding 
fireball. The time variability of the phenomena and the 
compact nature of the source suggest that the expand- 
ing wind has a bulk Lorentz factor of a few hundreds, 
a condition in principle sufficient to accelerate charged 
particles up to 100 EeV. In a more recent analysis taking 
into account the cosmological nature of the GRB distri- 
bution similar conclusion has been drawn, placing GRB's 
as one of the prominent sites of cosmic ray acceleration. 
Note that in such a scenario, UHECR sources are not 
visible since the detected cosmic rays come from various 
bursts and reach the Earth l ong after (10'^ t o 10^ years) 
the gamma ray burst itself ( (|Waxmanl . boOGl )). 

Direct observation of radio galaxies gives us their main 
characteristics in term of their radio luminosity (10"^^ — 
10^* ergs/sec), their size (10'^ — 10^ parsecs), their bright- 
ness morphology and the polarization level of the radio 
emission. From these parameters one can indirectly infer 
their mean magnetic field (of the order of 10 — 10'^ /uG) 
and kinetic power (10^^ — 10*^ ergs/sec). However the 
exact characteristics of the jets and in particular their 
Lorentz factor, density a nd composition are still under 
debate ( (|Massaghal l2008l) ). Among radio loud galaxies 
the Fanaroff-Riley radio galaxies of class II are of par- 
ticular interest because they combine a very powerful 
engine and relativistic blast wave (with Lorentz factor 
of the order 2 — 10) together with a relatively scarce 



environment. Hence, in the associated hot spots where 
the relativistic jet terminates, they not only satisfy the 
acceleration criterion but also the requirement that the 
accelerated particle does not lose all of its energy via 
r adiation or interactions on it s way out of the source 
fi Rachen and Biermann . 19931 )). Finally, invoking the 
sheared jet mechanisms where inductive acceleration can 
take place at the interface of the central spine and outer 
flow of the jet, accele ration of UHECR can take place in 
the jets themselves (( Lvutikov and OuvedL 2007)). 



B. Exotic top-down models 

One way to overcome the many problems related to 
the acceleration of UHECR is to introduce the existence 
of a new unstable or meta-stable super-massive particle. 
Its decay should produce quarks and leptons, which will 
result in a large cascade of energetic photons, neutrinos, 
and light leptons with a small fraction of protons and 
neutrons. In such a model no acceleration is required and 
cosmic rays are emerging directly with ultra high energy 
from the decay cascade. Hence their name of top-down 
models. 

For this scenario to produce observable particles above 
50 EeV three conditions must be met: 

• The decay must occur in recent time, i.e. at distances 
less than about 100 Mpc 

• The mass of this new particle must be well above the 
observed highest energy (100 EeV range), a hypothe- 
sis well satisfied by Grand Unification Theories (GUT) 
whose scale is around 10^ — 10^ EeV. 

• The ratio of the volume density of this particle to its 
decay time must be compatible with the observed flux. 
Two distinct mechanisms may produce such energy re- 
lease. 

• Radiation, interaction or collapse of Topological De- 
fects (TD), producing GUT particles that decay in- 
stantly. In those models the TD are leftovers from the 
GUT symmetry breaking phase transition in the very 
early universe. However very little is known about the 
phase transition itself and about the TD density that 
survives a possible inflationary phase, and quantitative 
predictions are usually quite difficult to rely on. 

• Super-massive metastable relic particles from some pri- 
mordial quantum field, produced after the inflationary 
stage of our Universe. Lifetime of those relics should be 
of the order of the age of the Universe and must be guar- 
anteed by some almost conserved protecting symmetry. 
It is worth noting that in some of those scenarios the relic 
particles may also act as non-thermal Dark Matter. 

In the case of TD the flux of UHECR is related to their 
number density and their radiation, collapse or interac- 
tion rate, while in the case of massive relics the flux is 
driven by the ratio of the density of the relics over their 
lifetime. 
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The very wide variety of topological defect mod- 
els together with their large number of parameters 
makes them difficult to review in detail. Many au- 
thors have addressed this field. Among them let us 



ment i on ([Vilenkin and Shellardl . Il995f) and (jVachaspati , 
19971 Il998r ) for a review on TD formation and inter- 
action. For a review o n experimental signatu r es see 
(iBhattacharied. Il998h . (iBhattachariee and Sigl I2OOOI ) 
and (jBerezinskv et al . 19981 ) ' 

Basic principles ruling the formation of TD in the early 
universe derive from the current picture on the evolution 
of the Universe. Several symmetry breaking phase transi- 
tions such as GUT ^ H ... =^ SU{3) x SU{2) x [/(I) 
occurred during the cooling. For those "spontaneous" 
symmetry breakings to occur, some scalar field (similar 
to the Higgs field generating masses to elementary par- 
ticles) must acquire a non vanishing expectation value 
in the new vacuum (ground) state. Quanta associated 
to those fields have energies of the order of the symme- 
try breaking scale, e.g. 10^^ — 10^^ GeV for the Grand 
Unification scale. 

During the phase transition process non causal regions 
may evolve towards different states in such a way that at 
the different domain borders, the Higgs field keeps a null 
expectation value. Energy is then trapped in a TD whose 
properties depend on the topology of the manifold where 
the Higgs potential reaches its minimum. Possible TDs 
are classified according to their dimensions: magnetic 
monopoles (0-dimensional, point-like); cosmic strings (1- 
dimensional) ; a sub- variety of the previous which car- 
ries current and is supra-conducting; domain walls (2- 
dimensional) ; textures (3-dimensional). Among those, 
only monopoles and cosmic strings are of interest as pos- 
sible UHECR sources. 



Supermassive 
possible source 



relic particles may be 
of UHECR ( (jBerezinskv 



a nothe r 
19991) . 



(jBhattachariee and Their mass should 



be larger than 10^^ GeV and their lifetime of the order 
of the age of the Universe since these relics must decay 
now (close by) in order to explain the UHECR flux. 
Unlike strings and monopoles, relics aggregate under the 
effect of gravity like ordinary matter and act as a (non 
thermal) cold dark matter component. The distribution 
of such relics should consequently be biased towards 
galaxies and galaxy clusters. 

Regardless of the details and dynamic of the topologi- 
cal collapse or of the massive particle decay the cascade 
that is produced will contain, possibly among many other 
things, quarks, gluons and leptons. Those particles will 
in turn produce far more photons and neutrinos than any 
type of nucleons. Hence, in all conceivable top-down sce- 
narios, photons and neutrinos dominate at the end of the 
hadronic cascade. This is the important distinction from 
the conventional acceleration mechanisms. The spectra 
of photons and neutrinos can be derived from the charged 



and neutral pion densities in the jets as: 
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where Ejct is the total energy of the jet (or equivalently 
the initial parton energy). Since 'i>^±(e,t) ~ 2$^o(e,t), 
photons and neutrinos should have very similar spectra. 
These injection spectra must then be convoluted with the 
transport phenomena to obtain the corresponding flux 
on Earth. In particular the photon transport equation 
strongly depends on its energy and on the poorly known 
univers al radio backg round and extragalactic magnetic 
fields f ijStanevl . 120101' )'). Nevertheless, if top-down sce- 
nario dominates the UHECR production above a certain 
energy the photon fraction should become very large. 
However, the recent observations, in particular of the 
Auger observatory, showed the contrary. This has con- 
siderably reduced the possibility for such models to be 
the source of UHECR. 



V. UHECR DETECTORS 

A. Older experiments - AGASA 

The AGASA (Akeno Giant Air-Shower Array) is the 
largest air shower array of the previous generation of de- 
tectors. AGASA covered 100 km^. It consisted of 111 
scintillator counters of area 2.2 ni^ at an average distance 
of 1 km from each other. Initially AGASA was divided in 



four branches that operated individually ( (jChiba et al. 



1992|)). In 1995 the d ata acquisition system was improved 
( (jOhoka et al\ . \l99l\\ ) and the four branches were unified 
in a single detector. This increased the effective detector 
area by a factor of 1.7 to reach 100 km^. 

Each AGASA station was viewed by a 125mm pho- 
tomultiplier and had a detector control unit (DCU) that 
controled the high voltage of the PMT, adjusted the gain 
and recorded the timing and pulse hight of every signal. 
The stations were connected by two optical fibers. One 
of them was used to send commands to the DCU. The 
other reported to the center triggers, shower data, and 
monitor data. AGASA operated for more than 12 years. 

Muon detectors of sizes from 2.4 m^ to 10 m^ were 
installed at 27 of the detector stations. In the Southeast 
corner of AGASA was the Akeno 1 km^ array which has 
been in operation since 1979. It was a densely packed 
array with detector separations from 3 to 120 m. Akeno 
has studied the cos mic ray energy spect rum from 3x 10^** 
eV to 3x lO^* eV ( IfNaeano et all Il984l) '). 
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B. HiRes 

The Hires observatory was a much improved foUow-up 
of the pion eer and very successful fluorescence detector 
Fly's Eye f (|Baltrusaitis et g/J . [l985h V 

Also constructed by the University of Utah, the ob- 
servatory was comprised of two air fluorescence detector 
sites separated by 12.6 km f (|Abbasi et gZ.l . l2004Ll2005c[ )). 
It was located at the U.S. Army Dugway Proving Ground 
in the state of Utah at 40.00°N, 113°W, at atmospheric 
depth of 870 g/cm^. The two detectors, referred to as 
HiRes-I and Hires-II operated on clear moonless nights 
with a effective duty cycle for physics data of about 10% 
typical for fluoresce nce detectors. 

The HiRes-I site (( Abu-Zavvad et al . 1999() ) consisted 
of 21 telescope units, each equipped with a 5 m^ spher- 
ical mirror and 256 phototube pixels at its focal plane. 
Each telescope covered an elevation range of 14° between 
3° and 17° and 360° in azimuth. The phototubes were 
equipped with sample and hold electronics which inte- 
grated the fluorescence signal within a 5.6 fis window. 
This was enough to contain the shower signal but also 
because of the limited elevation range of the detectors 
did not allow to extract the shower geometry from Eq. 1141 
alone. HiRes-I was in operation from June 1997 up until 
Aprfl 2006. 



12.6 km 



HiRes 1 




HiRes 2 



FIG. 12 Sketch of the HiRes fluorescent experiment. Each 
rectangle represents a fluorescence telescope including a mir- 
ror and a camera. Each site of the HiRes detectors has a 
nearly full azimuth coverage and site 2 which consists of two 
rings of mirrors covers elevation from 3° to 30° while site one, 
with a single ring, covers elevation from 3° to 16°. 

The HiRes-II site was completed at the end of 1999. 
Detectors were similar to those of HiRes-I but with twice 
as many mirrors organized in two rings covering eleva- 
tion from 3° to 31° and still 360° in azimuth. Moreover 
the HiRes-II phototubes were equipped with fast FADC 
electronics which sampled the shower signal every 100 ns. 
This allowed the reconstruction of the shower geometry 
from tim ing alone (Eq . I14II w ith a precision of about 5 



degrees f ()Abbasi et g/.l . l200' 



Although the two detectors of HiRes could trigger and 
reconstruct events independently, HiRes was designed to 
measure the fluorescence light stereoscopically. Stereo- 
scopic mode allows the reconstruction of the shower ge- 
ometry with a precision of 0.4° and provides very valuable 
information and cross checks about the atmospheric con- 



ditions at the time of the event. HiRes-I and II took data 
until April 2006 for an accu mulated exposure in stereo- 



scopic mode of 3460 hours C dAbbasi etaH . 120091) '). On 
the other hand, the "monocular" mode had better sta- 
tistical power and covered a much wider energy range. 

In monocular mode, the geometry of the HiRes-I 
events was calculated using an expected form of the 
shower development in addition to Eq. [14] (the PCF 
technique) . The shower profile was assumed to be de- 
scribed by the Gaisser-Hillas parameterization which is in 
good agreem ent with HiRes measurements and detailed 



simulations ((lAbu-Zayvad et all 120011 : iKalmvkov et al 
19971 : ISong et allhOOd) ). Significant contamination froi 



the forward-beamed direct Cherenkov light degraded its 
reliability and tracks with xo > 120° or with large 
Cherenkov fraction (as estimated from Monte Carlo sim- 
ulation) were rejected. Monte Carlo studies showed that 
the RMS energy resolution for this method was better 
than 20% only at the highest energies (above 10^^'^ eV). 

For monocular reconstruction, from either HiRes-I or 
Hires-II, the aperture is energy and composition depen- 
dent and must be evaluated by Monte-Carlo simula- 
tions. The HiRes collaboration made extensive and de- 
tailed simulation of both the atmospheric cascade and 
their detector and studied the systematics uncertainty in 



the e stimation of the monocular aperture ( ([Abbasi et al. 



20071) ') .The stereo aperture was determined by the re- 



quirement that the Monte Carlo events trigger both tele- 
scopes. Because of the better reconstruction of stereo 
events the quality cuts for them were not as strict and 
the stereo aperture is somewha t higher above energy of 



1019-7 gv ( (jAbbasi et a/.l . l2009l )). It is smaller for events 



below 10i«-5 eV. 

To determine the correct shower energies, the air fluo- 
rescence technique requires accurate measurement and 
monitoring of the absolute gain of the telescope. In 
HiRes, two methods of calibration were used. One pro- 
vided nightly relative calibration and used a YAG laser 
connected to two mirrors, the other relied on a stable and 
standard light source and provided monthly absolute cal- 
ibration. The pulses from a YAG laser were distributed 
to 2 mirrors via optical fibers. They provided a nightly 
relative calibration. Relative phototube gains were sta- 
ble to wi thin 3.5% and the absolute gains were known 
to 10% ( (|Abu-Zavvad et g/.l . lioOOal) ). Fluorescence light 
from air showers is also attenuated by molecular diffusion 
(Rayleigh) and aerosol scattering. While the former is 
approximately constant, the aerosol concentration varies 
rapidly with time. At HiRes (likewise at the Auger obser- 
vatory) the aerosol content was measured by observing 
scattered light from steerable laser systems. 

The fluorescenc e yield has been mea sured 

by (jKakimoto et aLl . ll996[ ). by ([Nagano et all \2004. and 
more recently by the AIRFLY collaboration ( (|Ave et ali . 
2008 ))_. A review of those measurements is available 



in ( Argueros et ali 20081) . The HiRes collaboration used 
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the f luorescence spectrum compiled b y (IBunner et al. 



196?') and normalized it to the yield of ( Kakimoto et al. 
1996). 

For both HiRes-I and HiRes-II events, the photo- 
electron count was converted to a shower size at each 
atmospheric depth, using the known geometry of the 
shower, and correcting for atmospheric attenuation. The 
reconstructed profile was integrated over the atmospheric 
depth. The integral was then multiplied by the average 
energy loss per particle to give the visible shower energy. 
A correction (about 10%) for the invisible energy, carried 
off by non-observable particles, was applied to give the 
total shower energy. 

The HiRes data contains two events at or above 
10-^° eV, measured at 1.0 and 1.5x10^° eV. Assuming a 
purely molecular atmosphere a lower energy limit of 0.9 
and 1.2x 10^° eV was obtained for these events. The flux 
values were on average 13% lower than the stereo spec- 



trum reported by Fly's Eye collaboration ([Bird et al. 



19931 ). This difference can be explained by a 7% offset in 
the energy calibration, well within the uncertainty of the 
two experiments. 




FIG. 13 The Pierre Auger observatory at the end of March 
2009. Individual white dots represent Cherenkov tanks, while 
gray ones are un-equiped positions. A denser (infill) area is 
visible in the upper left. Big white dots at the periphery 
of the array are fluorescence detector sites with the field of 
view of individual telescope given by the radial white line. 
Also shown is the Central Laser Facility (CLF) used for FD 
calibration and atmospheric monitoring purpose. 



C. Auger 

The Pierre Auger Observatory is the largest operating 
cosmic ray observatory ever build. It is based on the 
hybrid concept where both fluorescence and surface array 
detection techniques are used and combined. 

The Southern site of the Auger observatory is located 
in the "Pampa Amarilla" region (35.1°-35.5° S, 69.0°- 
69.6° W and 1300 -1400 m a.s.l.) of the p rovince of Men- 
doza, Argentina (( Abraham et ali\2004 )). Construction 



was completed in 2008 but stable data taking started as 
early as the beginning of 2004 when Auger already had 
100 detectors, covering and area in excess of 150 km^, 
installed in the held . The arrangement of the detectors 
is shown in Fig. 1131 



1. Surface Array 

The surface array (SD) of Auger South is composed 
of 1600 water Cherenkov tanks, distributed on a trian- 
gular grid of 1500 m. It covers a total surface area of 
3000 km^. Each tank is equipped with three photomul- 
tiplier tubes to measure the Cherenkov light, a data ac- 
quisition (DAQ) and front-end electronic (EE) card for 
control and trigger, a solar panel and two batteries for 
power, a GPS receiver for the time tagging, and a cus- 
t om radio emitter /rece iver for trigger and data transfer 



(( Allekotte et al . 2008f )). A central site, located on the 



Southwest corner of the array hosts the central DAQ sys- 
tem (CDAS), including the central trigger processors and 
the permanent data storage area. 



The SD has a 100% duty cycle, and a well defined 
purely geometrical aperture (cx cos 9) above trigger sat- 
uration at 3x10""^® eV. The coverage is largely uniform in 
right ascension. Modulation in the event rate due to the 
atmospheric conditions are at the level of 2% for daily 
modulation and about 10% for seasonal ones. Those ef- 
f ects have been carefully studied and can be corrected for 
r (|Abraham et a/.l . l2009a|) V 

The water tanks are 1.2 m in height and are mainly 
sensitive to muons, electrons, positrons and photons. A 
vertical GeV muon hitting the tank deposits an energy 
of about 240 MeV, to be compared to few tens of MeV 
for an average electron. The unit for the shower signal 
is a vertical equivalent muon (VEM). This allows for an 
in situ calibration of the PMT gain based on the rate 
of atmospheric muons. The gain is adjusted so that a 
single ADC count corresponds to about 1.5 MeV. Local 
triggers are adjusted to a rate of about 20 Hz for a simple 
threshold trigger and a few Hz for a more sophisticated 
time over threshold (counting time bin s over a certain 



threshold within a given time window) (( Abraham et al. 



2010eh ). Local triggers are sent to CDAS where space- 
time coincidences of at least three tanks are required to 
trigger the upload and permanent storage of the full event 
data. 

The large sensitivity to muons and the height of the 
individual tanks allows the Auger array to have excellent 
sensitivity to horizontal showers be them from hadronic 
origin or from neutrinos. 

The shower signal is sampled at a rate of 40 MHz. The 
analysis of its time structure allows for example to iden- 
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tify the presence of an electromagnetic component in the 
ground signal, at appropriate distances from the core to 
identify the short high pulse from the individual muons 
and to count them and to calculate signal shape param- 
eters such as the rise time. Those information allow to 
efficiently distinguish neutrinos from the hadronic back- 
ground in nearly horizontal showers and photons. Addi- 
tionally they also allow to construct hadronic mass sen- 
sitive parameters. Timing information is obtained from 
GPS receiver functioning in position hold mode. The 
absolute time resolution is about 10 ns, combined with 
the sampling of the shower front and of the FE electron- 
ics which allows for an angular reso l ution of bet ter than 
1° above lO^^ eV f ljBonifazi et al\ . l2008l l2009h 'l. The 
aperture of the Southern Auger Observatory is energy 
independent when the surface array triggers and is deter- 
mined by the area of the SD and the maximum shower 
zenith angle (60°)used in the analysis. 

The lateral distribution functio n of the Auger tank 
signals is fitted to a NKG f ljGreisenl . mm and 
( Kamata and Nishimura . 1958f) ) function: 



f{r) 
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1000 m + 700 m 



/?+7 



(18) 

where S'looo is the adjusted normalization and where the 
exponent /3 is adjusted to the data using a second order 
polynomial in sec 9 whose coefficient is a linear function 
of ^1000 in VEM(e.g. a = ao+oi logio(5iooo/VEM). The 
exponent 7 is very close to zero. 



2. Fluorescence detector 

The fluorescence detector of Auger South (FD) is com- 
posed of 24 telescopes distributed in four sites installed 
at the periphery of the surface array and looking inward. 
Each telescope has a field of view of 30°x30° in elevation 
and azimuth. A set of six telescopes in each site covers 
180° in azimuth and observes the atmosphere above the 
ground array. This geometrical arrangement ensures full 
detection efficiency for showers in excess of 10^ ^ eV over 



the entire surface of the array ( (lAbraham et aLl . l2010a|) ). 

In each telescope the optical system is composed of an 
entrance filter selecting the UV light, an aperture and 
corrector ring maintaining a large aperture while reduc- 
ing spherical and eliminating coma aberrations, and a 
3.6 m diameter mirror illuminating a camera composed 
of 440 PMT tubes. Each tube has a field of view of 
1.5°xl.5°. 

Triggering is done at the hardware level of each cam- 
era for the first (pixel) and second (alignment) levels. 
A third level trigger (TLT) is implemented in software 
mainly to reject lightning events and random alignments. 
Each TLT is then processed at the fluorescence site level 
to merge all the telescope information and to send via 



CD AS, a preliminary shower direction and ground im- 
pact time to the surface array. The information from the 
tanks closest to the shower core is retrieved for showers 
that do not independently trigger three tanks. Together 
with additional fiducial cuts this hybrid trigger is fully ef- 
ficient above 10^^ eV. Above this energy, the FD trigger 
is always accompanied by at least one station, indepen- 
dent of t he mass and direction o f the incoming primary 
particle (( Abraham et al . 2010cil )). 

Event reconstruction proceeds in two steps. First the 
shower geometry is found by combining information from 
the shower image and timing measured with the FD with 
the trigger time of the surface dete ctor station that has 
the largest sig nal ( (iMostafal . l2007h 'l. From this timing 
information it is possible to break the degeneracy in- 
trinsic to equation 1141 Therefore the hybrid approach 
to shower observation enables the shower geometry and 
consequently the energy of the primary particle to be 
determined accurately. The Auger collaboration uses a 
fiuorescence yield in air at 293 K and 1013 h Pa from 
the 337 nm band of 5.05±0.71 photons/MeV of energy 
deposited take n from the me a surem ents of Nagano and 
collaborators (( Nagano et al . 20041 )). The wavelength 



and pressure dependence of the yield adopted follows the 



meas urements of the AIRFLY collaboration ( (|Ave et al. 
20081 )). Note that th e fiuorescence yi el d use d by the 



HiRes collaboration (^( Kakimoto et al\ 19961)) in the 



same conditions is 5.4 photons/MeV ( ([Argueros et 



20091) 1 so very close to the value used by Auger. 

In the second step light attenuation from the shower 
to the telescope is estima ted and all contrib uting light 
sources are disentangled (( Unger et al . 20081 )). A great 



deal of effort is spent by the Auger collaboration to accu- 
rately monitor the atmospheric transparency and main- 
tain the absolute calibration of the telescopes. An ex- 
tensive set of instruments is installed and operated at 



the A uger site for this sole purpose ( ([Abraham et 



2010bl )). Finally the profile of energy deposition of the 
shower is reconstructed using a Gaisser-Hillas functional 



form ( ([Abraham et a/.l . lioiOal) ) 



The reconstruction accuracy of hybrid events is much 
better than what can be achieved using SD or FD data 
independently. For example, the angular and energy res- 
olution of hybrid measurements at 1 EeV is better than 
0.5° and 6% respectively compared with about 2.5° and 
20% for the surface detector alone. 



VI. COSMIC RAY ENERGY SPECTRUM 
A. The end of the cosmic ray spectrum 

In 2008 the HiR e s Coll aboration published a pa- 



per ( (jAbbasi et al\ . l2008a[ )) with a title emphasiz 



ing the experimental proof of the GZK suppression. 
Soon after that the Auger Collaboration confirmed the 
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o bservation of the end of the co smic ray spec trum 



19661) and 



f (|Abraham et all l2008bD. In 1966 (iGreisep . . _ _ 
independently ( Zatsepin and Kuzminl . 1966 ) predicted 
that the cosmic ray spectrum wih end at several times 
10^® eV because of the interactions of the UHECR with 
the microwave background (CMB). Although the energy 
of the CMB photons is very low, the center of mass en- 
ergy of these interactions is enough to produce pious and 
cause high energy loss for these particles that decreases 
their flux. 

The importance of these observation is that the pre- 
viously biggest air shower array, AGASA, has observed 
11 events abov e 10^° eV and no de crease above the pre- 
dicted cutoff (( Takeda et al . 1998h ). A re-evaluation of 



the energy assignment of AGASA based on a la rger data 
set was published later in (i Takeda et all . l2003l ). The en- 
ergy determination of AGASA was tied up to the particle 
signal at 600 meters from the shower axis S'o(600). The 
Monte Carlo calculations suggested that the primary en- 
ergy is 

E = 2.17 X 10i °^£;o(600) eV. 

Since the detectors of AGASA were scintillator coun- 
ters the signal is produced by the shower electromag- 
netic component with a small contribution of the shower 
muons. The ultrahigh energy events published by 
AGASA provided the inspiration for the exotic 'top- 
down' models of these particles. 

The HiRes energy spectrum was based on monocular 
observations with the two fluorescent telescopes HiRcs-I 
and II that were operated respectively from 1997 to 2005 
and from 1999 to 2004. Some corrections were made 



on th e previous HiRes spectrum release ( (jAbbasi et i 
20041) 1. A special attention was paid to the detector 



calibration and the atmospheric conditions which were 
studied by standard meteorological methods and by the 
observation and analysis of laser shots from different lo- 
cations surrounding the two detectors. Most of the high- 
est energy events were observed with HiRes-I. The exact 
shape of the spectrum depends strongly on the calcula- 
tion of the aperture for the two telescopes. The HiRes 
found two breaks in the cosmic ray spectrum - one at 
energy lois-^^ito.os (^.-^^ cosmic ray ankle) and another 
at I0i9-75±0 04 ^j^g Q2K cutoff. The spectral index 
between the two breaks is 2.81±0.03 and after the cut- 
off is 5.1 ±0.7. The spectrum is consi stent with various 



models and in particular the model of (jBerezinskv et al 



[2005) with pure proton composition. An important pa- 
rameter is Ei/2 where the cosmic ray flux is one half of 
what it should h ave been without the GZK effect. i?i/2 
was predicted by (jBerezinskv and Grigor'evaill988l ) to be 
lO^^ '^s and HiRes measured I0i9-73±o.o7^ rpj^^ statistical 
significance of the cutoff is more than 5a. 

The energy spectruni derive d by the Auger collabo- 
ration ( ([Abraham et al . 2008b[) ) is shown with those of 
HiRes and AGASA in Fig. [TH As well as HiRes, Auger 



> 

0) 



(0 
'(0 

CM 

'e, 

z 

■a 

CO 

LU 



10'- 



AGASA 
HiRes 1 
HiRes 2 
Auger 



♦ ♦ ♦ 



n 



10" 



10'^ 

Energy (eV) 



10" 



FIG. 14 Cosmic ray spectrum as presented by AGASA, HiRes 
and Auger. 



observes the shower profile with its fluorescent detectors. 
They, however, have a live time of 13% compared to that 
of the surface detector. Taking full advantage of its hy- 
brid design the Auger collaboration decided to correlate 
the energy derived from the fluorescence observations to 
the shower signal at 1,000 meters from the shower core 
(S 1000) 5 which is the least sensitive to the cosmic ray 
composition. To account for the angular dependence of 
this quantity it was corrected to the median angle of 38° 

c38 



n using the constant intensity curve f (|Hersil et al 
19611 )1 observed by the surface array. The constant in- 



tensity curve is a study of the change with angle of the 
signal threshold above which cosmic rays arrive with con- 
stant rate, i.e. study of the shower absorption in the 
atmosphere. 

The correlation between FD and Siooo was studied in 
high quality hybrid events, that were seen in both the 
surface and the fluorescent detectors. The correlation 
showed that 



Efd = 1.49±0.06±0.12(syst)xS'- 



1.08±0.01±0.04 
1000 



xlO^^ eV 



(19) 

The Auger collaboration then used the surface detector 
statistics to produce the energy spectrum. The uncer- 
tainty in the energy reconstruction by the fluorescent 
telescopes was estimated to 22% and the width of the 
observed correlations was consistent with the statistical 
uncertainty of both measurements. The surface detec- 
tor exposure for this publication was twice that of HiRes 
and four times higher than AGASA. Since the fiuorescent 
energy measurement does not depend on the hadronic in- 
teraction model used in the analysis, such an estimation 
of the spectrum was considered to be model independent. 

The Auger spectrum has a slightly different shape 
in addition to the energy assignment. From 4x10^* 
eV to 4x10^^ eV the slope of the spectrum is 
2.69±0.02±0.06(syst) and above it is 4.2±0.4±0.06(syst). 
A single power law for the whole data set is rejected at 
6cr level. 
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The measurements of the cosmic ray s pectrum were 
exten ded by using stereo events in HiRe s ((lAbbasi et al. 



20091)') and hybrid events in Auger ( () Abraham et i 
2010d|))). Stereo events are reconstructed much more 



precisely than monocular ones and they confirmed the 
previously measured spectrum. Auger used hybrid events 
to extend the spectrum to lower energy. The Auger ex- 
posure at the time of this last publication was 12,790 
km^yr.sr. All measured spectra from HiRes and Auger 
are shown in Fig. 1151 
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FIG. 15 Cosmic ray spectrum as measured by HiRes and 
Auger. 

The new Auger spectrum is a bit flatter than the older 
one with an index of 2.59±0.02 between the breaks and 
4.3±0.2 above that. i;i/2 value is lOi^-^iio os. The differ- 
ences in the interpretation of the HiRes and Auger spec- 
tra are significant. The Auger spectrum can be explained 
by several different models some of which include mixed 
chemical composition at acceleration in the sources. The 
end of the cosmic ray spectrum measured by Auger is 
consistent with the GZK effect. 



B. Cosmic ray energy loss in propagation 

In addition to the adiabatic energy loss because of the 
expansion of the Universe, there are two important en- 
ergy loss processes for protons: pion photo-production 
interactions and e^e~ pair production (BH) interactions 
identical to the pair production interactions of 7-rays in 
the nuclear field. The average interaction length \ph for 
interactions with the CMB is the inverse of the product 
of the interaction cross section a^h and the photon den- 
sity n. For (jph = 10~^* cm^ and n = 400 cm~'^, \h = 
8.3 Mpc. 

Heavy nuclei lose energy in photo-disintegration (spal- 
lation) processes when the center of mass energy exceeds 
the giant dipole resonance. Since less energy is required 
in the center of mass, the cross section is higher but the 
energy loss depends on the mass of the nucleus that loses 
one or two nucleons. The photo-production energy loss 



follows the same energy dependence as for protons but 
in the Lorentz factor space, i.e. in E/A units. The pair 
production cross section is a quadratic function of the 
charge of the nucleus Z . 

In the case of 7-rays the energy loss is due to the 77 — > 
e^e^ process. 

A photo-production interaction is possible when the 
center of mass energy of the interaction ^/s is higher than 
the sum of a proton mass rUp and a pion mass m^r. In 
the Lab system the square of the center of mass energy 
is 



2Epe{l 



(20) 



where e is the photon energy and 9 is the angle be- 
tween the proton and the photon. In a head on collision 
{cosO = — 1) with a photon of the average CMB energy 
(6.3x10^'* eV) the minimum proton energy is 



En 



4e 



(2mp -f m^r) 



10^0 eV. 



(21) 



There are many CMB photons with higher energy and 
the threshold proton energy is actually lower, about 
3x1019 eV. 

The cross section for pion photo-production was very 
well studied at accelerators. The highest cross section 
is at the mass of the A+ resonance (1232 MeV). At the 
peak of the resonance the cross section is about 500 ^b. 
The cross section decreases to about 100 /xb and then 
increases logarithmically. The neutron interaction cross 
section is very similar to the proton one. 

The CMB spectrum and density are also very well 
known, so the proton interaction length can be calcu- 
lated exactly. Since protons lose only a fraction of their 
energy {Kind ) , another quantity - the energy loss length 
Lioss — becomes important. The energy loss 

length is longer than the interaction length by 1/Kinei, 
by about a factor of 5 at threshold. At higher energy 
Kind grows and this factor is about 2. 

In the case of e~ ^e~ pair production 

19881) ') 



the addition of 



(jBerezinskv and Grigor'eva . 
two electron masses to the center of mass energy ^/s 
requires much lower proton energy and the process has 
lower threshold. The cross section is higher than aph, 
but the fractional energy loss is of order of me/ rap. The 
energy loss length has a minimum around 2x10^^ eV 
and is always longer than 1,000 Mpc. 

The last proton energy loss process is the redshift due 
to the expansion of the Universe. The energy loss length 
to redshift is the ratio of the velocity of light to the 
Hubble constant {c/Hq) and is 4,000 Mpc for i/o = 75 
km.s~iMpc~i. 

The energy loss length of protons in the CMB is shown 
in Fig. [HI 

The energy loss length for several nuclei i s also shown 
in Fig. [m as calculated by (jAllard et all l2005l ). The 
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minimum value of Lioss is significantly lower than that 
of protons but is achieved at higher energy: AxEp. Since 
only iron has similar L/oss to protons around 10^° eV it is 
considered the only nucleus that can compete with pro- 
tons in the chemical composition of UHECR. The effect 
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FIG. 16 Energy loss length for protons, nuclei and gamma 
rays. The heavier shading points at the proton and gamma- 
ray Lioas and the light one shows the contribuion of the e"'"e~ 
pair production. The adiabatic energy loss is not included. 

of propagation on the accelerated UHECR can not be 
calculated directly from the energy loss lengths shown in 
Fig. [16] because an accelerated nucleus changes its mass 
after the first photodisintegration. A code treating the 
propagation of nuclei should account for the energy loss 
of all nuclei and isotopes lighter than the injected nu- 
cleus. 

The process 77 — >■ e+e^ has a resonant character and 
the cross section peaks a.t E^e = 2nil where e is the 
ambient photon energy. For CMB this corresponds to 
Ej of 8x 10^* eV and the mean free path decreases with 
increasing Ej. For gamma rays of energy 10^" eV the 
relevant seed photon frequency is about 1 MHz - in the 
radio band. This creates some uncertainty in the esti- 
mates of the UHE 7-ray energy loss length because the 
density of the radio background at such frequencies is not 
well known. 

A different source of uncertainty in the 7-ray propaga- 
tion is the strength of the extragalactic magnetic fields. 
If they are negligible the electrons have inverse Compton 
interactions, whose interaction length is similar to that 
of the pair production, and generate a second generation 
of very high energy 7-rays. If, however, the magnetic 
fields are significant, electrons lose energy very fast on 
synchrotron radiation and the created 7-rays are in the 
MeV-GeV energy range. The energy loss distance on 
synchrotron radiation is 2.6E^^BZg Mpc, where Eis is 



the electron energy in units of 10^® eV and _B_g is the 
strength of magnetic field in nGauss. 



C. Formation of the cosmic ray energy spectrum in 
propagation 

Predictions of the shape of the cosmic ray spectrum 
requires much more than the energy loss in propagation. 
The necessary astrophysical input includes at least the 
following items: 

• UHECR source distribution 

• cosmic ray source emissivity 

• cosmic ray injection (acceleration) spectrum 

• maximum acceleration energy T^^rnax 

• cosmic ray chemical composition 

• cosmic ray source cosmological evolution 

that are not independent of each other. As an example 
we will discuss the formation of the proton spectrum in 
propagation. 
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FIG. 17 Contribution of different redshifts to the arrival spec- 
trum for E~'^ injection spectrum with no cosmological evolu- 
tion. The thick gray line shows the sum of the contributions 
from these five redshifts while the black line is result of a full 
integration. 

Figure [T71 shows the contribution to the observed UHE 
cosmic ray proton flux by sources located at different 
redshifts that inject protons on a E~^ spectrum with an 
exponential cutoff at 10^^'^ eV. One can see how the en- 
ergy loss increases the contribution to the 2-6x10"'^^ eV 
energy range after propagation to z=0.1. In models with 
cosmological evolution of the sources the effect is stronger 
and proportional to the strength of the source evolution. 

A simplification in such calculation is the assump- 
tion that sources are isotropically and homogeneously 
distributed in the Universe and the contribution of all 
sources are identical. In such a case the cosmic ray flux 
at Earth can be determined by an integration of the fluxes 
from different redshifts shown in Fig. [T71 In the case of 
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cosmological evolution of the sources the integral is 



A. Limits on the flux of neutrinos 



^° dt 

L{z)No{Eo)P{Eo,E', z)—dE'dz 



./n .IP, UZ 

(22) 

where L(z) is the cosmic ray source emissivity as a func- 
tion of redshift and No{Eo) reflects the injection spec- 
trum. P{Eq, E' , z) is the probability for a proton injected 
with energy Eq at redshift z to reach us with energy E' . 
The derivative dt/ dz depends on the cosmological model 
and is 



Possible showe r neutrino primaries may be the eas- 



dt 
dz 



1 



Ho{l + z 



and is simplified to (1 + z)-^/'^ / {Ho{l + z)) 
Einstein-deSitter Universe. 

It is important to note that the contribution of different 
redshifts depends not only on the cosmological evolution 
but also on the injection spectral index as the photopro- 
duction energy loss is a strong function of the injection 
energy. Since in steep injection spectra a larger fraction 
of the observed flux comes from lower primary energy 
(that do not change as much on propagation) the contri- 
bution of higher redshifts is larger. 

One can see in Fig.[T7]that even 2=0.05 does contribute 
to UHECR above GxlO^^ eV where the GZK cutoff is 
already present. Another explanation is that the cutoff 
is just the end of the acceleration p ower of the sources 
that does not much exceed 10^° eV fl Aloisio et al. . 20091 : 



Watson 



200?! )). The extragalactic magnetic fields can 
also be involved in the explanation. If they are high 
UHECR would scatter often and their real pathlength 
would be con siderably larger tha n the distance to the 
sources as in ( Stanev et al. . I2OOOI) . 



VII. CHEMICAL COMPOSITION OF UHECR 

One has to use the properties of the extensive air show- 
ers to identify the type of the primary particle whose 
interaction in the atmosphere has initiated the shower. 
Because of the high level of fluctuations in the shower 
development it is quite difficult to distinguish showers 
originating from different hadronic primaries on an event 
by event basis - it can only be done on a statistically 
significant set of showers. At lower energy, around the 
knee, the main parameter in composition studies is the 
ratio of the shower muon to electron components which 
increases with the primary nucleus mass. 

Another way is to study the shower longitudinal devel- 
opment. This is usually done by observing the depth of 
shower maximum Xmax with fluorescent detectors that 
can determine the atmospheric depth where the shower 
particles emit the highest amount of fluorescent light. 
The shower longitudinal development, as we will see in 
this section, can also be studied by the surface detectors. 



iest to identify fi Capelle et al . 1998f )). The reason is 



the many orders of magnitude difference between the 
hadronic and neutrino cross sections. If neutrinos in- 
teract in the atmosphere at all, they would interact very 
deep. It is more likely that they interact in the rock of 
the Earth. This was used by the Auger collaboration in 
order to set a limit on the flux of t neutrinos. Setting 
such a limit is equivalent to a limit on the total neutrino 
flux. Although z^T-'s are rarely produced in particle in- 
teractions, cosmic neutrinos oscillate in propagation to 
Earth. While at production the neutrino flavor ratio 
for the {^e'-v^'.Vr) is close to 1:2:0 after propagation it is close 



to 1:1:1. 

The tau neutrino detection idea fi Bertou et al. . 20021 ) 1 
is that in a small fraction of the solid angle, at zenith an- 
gles 6 between 90.1° and 95.9° tau neutrinos will graze 
the Earth, possibly interact and after escaping the earth 
the tau decay will generate a shower that can be seen 
by the shower array. The neutrino identiflcation is based 
on the different quality of vertical and almost horizontal 
showers. Vertical showers are young: they exhibit long, 
^ fjLS waveforms in the surface detectors. Old showers, 
after penetrating about two atmospheric depths, consist 
mostly of muons. The waveforms they generate in the 
surface detectors are much shorter, of order of 100 ns. 
If one detects an almost horizontal shower that has the 
waveforms of young showers it would mean that the pri- 
mary particle has interacted near to the detector and is 
most likely a neutrino. 

Atmospheric interaction of tau neutrinos are also espe- 
cially interesting_becanase^^ develop two show- 
ers (( Learned and Pakvasa . 1995f )). one when the r neu- 
trino interacts and the second one when the r lepton de- 
cays. Since the t energy loss is much lower than that of 
a muon, most of the neutrino energy is released through 
the T decay. The exact parameters for v^- shower identi- 
flcation are found with Monte Carlo calculations. 

The Auger collaboration is using two related parame- 
ters. The first one is the shower shape in the surface ar- 
ray. Because of the large zenith angle it obviously should 
be elongated. The collaboration defined shower length 
and shower width. The Monte Carlo calculations showed 
that there is no chance for a nucleus or a 7-ray to gener- 
ate a shower with length/width ratio higher than 5. 

The second parameter is the shower ground speed. If 
the shower is indeed horizontal it has to move with ve- 
locity equal to the speed of light. So they looked for 
showers with velocity between 0.29 m/rts and 0.31 m/ns. 
Only showers with RMS (ground speed) better than 0.08 
m/ns are included in the sample. These requirements, to- 
gether with the general requirement that the tank with 
the maximum VEM signal is surrounded by six active 
tanks significantly decrease the size of the sample. Any- 
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way, no such showers were found in the Auger statistics. 
The Auger exposure as a function of energy was deter- 
mined by Monte Carlo calculations. 

Using the statistics between January 2004 and Au- 
gust 2007 the Auger collaboration set an integral limit 
on the Ur flux between 2x10^'' eV and 2x10^^ eV of 
E^dN/dEr of 1.3x10-'^ GeV. This limit assumes that 
the U t energy spectrum is E ~^. In the same publica- 
tion (; ([Abraham et aLl . l2008dh ) the collaboration uses the 
exposure as a function of the neutrino energy to also 
give a differential limit. Extending the statistics to April 
2008 the Auger collaboration decreased the integral limit 
o f EldN/dEr of to ab out BxlQ-s GeV.cm-^s-isr"! 



( (| Abraham et aLl . l2009c[ )) for the same fiat v-r spectrum. 
This limit is shown with a gray line in Fig. 1251 The in- 
tegrated Auger limit is competitive with the limits set 



by th e neutrino telescope AMANDA ( (jAchterberg et at. 
20071) '! at lower energy. 



The HiRes Collaboration has also set limits on the 



fluxes of tau and e l ectron neutrinos ( (|Abbasi et at. 



2mm . iMartens et all 120071 )). The first limit is based 
on the same assumptions as the Auger one. The HiRes 
collaboration has simulated r-neutrino induced showers 
hitting the Earth with elevations between 10° and -10° 
with an account of the topography of the detector. Af- 
ter detection simulation they obtained 6,699 monocular 
triggers and 870 stereo ones. Then the collaboration an- 
alyzed with some quality cuts simulated and real data 
events events with reconstructed zenith angles between 
88.8° and 95.5°. The data sample yielded a total of 134 
events that happened to be laser events, which passed 
the cuts because of the light scattering near the ground. 
Thus they were left with no neutrino candidates. 

For the limit on electron neutrino events only the 
HiRes 2 detector data were used because of its superior 
reconstruction. HiRes looked at upward going showers 
with zenith angles above 105°. Lower zenith angles do 
not yield more events because of the neutrino absorp- 
tion in the Earth. The basis of the search is the fact that 
high energy electrons have much lower energy loss at high 
energy and es pecially in dense materials becau s e of the 
LPM effect ( (jLandau and PomeranchuH 119531 : iMigdai 
19561) ). After electron neutrino charge current interac- 
tions in the ground the generated electrons would not lose 
much energy and may produce upward going air show- 
ers. The simulations were done for neutrinos of energy 
exceeding 10^^ eV that generate horizontal and upward 
going air showers with more than 10*" particles at maxi- 
mum. These showers were then treated with the HiRes 
detector simulations and compared to experimental data. 
No neutrino candidates were observed. The combination 
of the two searches reduced the integral neutrino flux 
limits in the three decades above 10^^ eV to 3.8xl0~^, 
9.7x10-6, and 4.7x10"^ GeV cm-^sr-^s-^ 



B. Limits on the fraction of gamma-rays 

The limit on the fluxes of ultrahigh energy neutrinos 
is astrophysically important because it is related to the 
dynamics of the systems where UHE cosmic rays are ac- 
celerated and the importance of hadronic interactions in 
such objects. The limit on the fraction of photons in 
UHECR determines the general origin of the highest en- 
ergy particles in the Universe. All top-down models of 
the UHECR origin predict that 90% of these particles 
are 7-rays and neutrinos. There have been previous lim- 
its from several UHECR air shower arrays for energies 
above 10^® eV but none of them was less than 10%. 

Auger has the advantage of being a hybrid array. A 
fluorescent detector trigger with a single surface detector 
trigger improves significantly the shower reconstruction 
and lowers the detection threshold. The Auger collab- 
oration set limits on the fraction of 7-rays at energies 
a bove 2, 3, 5, and lOxlO^s eV (2, 3, 5, and 10 EeV) 
( (| Abraham et all l2007al l2009el )). The limit is based on 
the measurement of the shower depth of maximum Xmax- 
Because of the low secondary multiplicity in electromag- 
netic interactions UHE 7-ray induced showers reach max- 
imum much later than proton showers. In addition at 
energies above 10^^ eV the LPM (Landau-Pomerntchuk- 
Migdal) effect, which significantly decreases the pair- 
production cross section, starts becoming important and 
increases Xmax even more. The key in such a measure- 
ment is to make certain that the event selection is not 
biased versus showers with deep X^ax- 

There are many cuts that are applied to the detected 
showers that exclude the possible biases. Auger requires 
that the tank with the largest signal is less than 1.5 km 
from the reconstructed shower axis and the time differ- 
ence between the fluorescent and tank signals is small. 
Another requirement is that the shower Xmax is observed 
in the telescope field of view. The minimum angle be- 
tween a fluorescent pixel and the shower direction has to 
be above 10° to exclude Cherenkov light contamination. 
The shower zenith angle has to be higher than 35° and 
the distance of the shower core to the telescope less than 
24 km. 

All these different cuts, together with the requirement 
for a good reconstruction in the fluorescent detector, de- 
crease the total statistics above 2 EeV to 2063 events. 
Eight of these events have Xmax consistent with possible 
photon showers. Using the 95% confldence value for the 
number of photon candidate events and the systematic 
uncertainty of 22% in the energy estimate and 11% in 
Xmax the Auger collaboration arrives at an 95% upper 
limit on the fraction of photon showers above 2 EeV of 
3.8%. At higher energy bins the number of photon can- 
didate showers is 1, 0, and 0, but the statistics are lower, 
1021, 436, and 131 events respectively. This leads to 
95% proton fraction estimates of 2.4%, 3.5%, and 11.7% 
as shown in Fig. [TH] 
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FIG. 18 Fraction of photon showers above certain energy de- 
termined by the hybrid analysis (labeled h) and the surface 
detector (sd). The shaded area shows limits set by previous 
experiments. 



At higher energy the photon fract ions are calculated 
using the surface detec tor results ((IRisse and Homola 
20071 )). This analysis ( (| Abraham et al 



2008d)) is very 



that were very different from those of the detected air 
showers from 2004 to the end of 2006. The data set used 
includes 2761, 1329, and 372 showers above 10, 20, and 40 
EeV. There were no 7-ray candidates in either bin, while 
570, 145, and 21 showers are certainly nuclear showers. 
Since zero events per bin corresponds to less than 3 events 
at 95% confidence level, the fraction of photon showers 
was calculated to 2.0%, 5.1% and 31% in the three bins. 

These results demonstrate that the top-down models 
are not responsible for the production of the majority 
of UHECR. There is still a little space remaining at the 
highest energies but in 5 years of operation Auger South 
will be able to bring down these limits if the current trend 
continues. 

It is indeed true that the differences between different 
hadronic interaction models do not play much of a role in 
the photon showers Monte Carlo and these analyses are 
mostly independent of the details of those interactions. 



interesting because it uses surface detector shower char- C. Depth of maximum data and their interpretation 



The idea of such analysis was first developed for the 
i nclined showers detected by the Haverah Park array 
( (|Ave et all l2002l )l When the air shower particles hit 



the surface array detectors the shower front is not flat, 
it has a certain curvature, i.e., the shower particles away 
from the shower axis arrive later than those close to it. If 
the shower curvature is assumed to be spherical (which is 
an oversimplification) the delay is proportional to / H 
where H is the altitude of the particle production and r 
is the distance from the shower axis. This means that in 
early developing showers the shower curvature is smaller 
than in later developing ones. The radius of shower cur- 
vature Rc is the first parameter that can be measured by 
the surface array. 

The second parameter is the width of the shower front, 
i.e. the time that it takes the shower particles to arrive at 
the surface array. The spread of the arrival time at cer- 
tain distance r from the core also increases with the depth 
of Xmax- This could be measured by the surface array 
as the arrival time at a fixed distance from the shower 
core T]^/2, which is defined by Auger as the time in which 
10% to 50% of the signal arrives at a detector. Since 
photon showers develop deeper in the atmosphere than 
nuclear showers, one can identify them by their large Rc 
and Ti/2- For photon showers one can use Monte Carlo 
calculations that have the advantage to depend very lit- 
tle on the hadronic interaction models used in nuclear 
shower calculations. The only remaining problem is the 
energy assignment of the photon showers - the fluores- 
cent detector value is used for hadronic showers. Auger 
developed an estimate that gave them 25% accuracy for 
photon showers. 

The Monte Carlo calculations of photon showers gener- 
ated values of these two parameters for all zenith angles 



The application of Heitler's toy shower model to the 
shower longitudinal development demonstrates its depen- 
dence on the mass of the primary particle. With the 
advent of fluorescence detectors the measurement of the 
shower depth of maximum, Xmax, quickly became a ma- 
jor component of the cosmic rays composition studies. 
An important parameter is the shower elongation rate 
Dio, the relation of which to the c hanges in the cosmic 
ray c omposition is discussed first in ( Linslev and Watsonl . 
19811) . 

The first analysis of the X^ax energy dependence with 
fluorescent dete ctor data was done with the Fly's Eye. 

which used only 



This analysis ( (jGaisser et al. 



1993D), 

two chemical components - H and Fe, showed a trend of 
increasing the proton fraction in cosmic rays of energy 
above 10^^ eV. 

In 2005 the HiRes Collaboration published an anal- 
ysis of t he UHECR composit ion from Xmax measure- 



ments (( Abbasi et al. . 2005bl )V The data sample in- 



cluded 553 events of energy above 10^* detected in stereo 
by both fluorescent detectors during 20 months from 1999 
to 2001. The sample is relatively small because of the dif- 
ferent cuts made on the total event sample. The first set 
of cuts are related to the atmospheric conditions. About 
3/4 of the total sample had hourly data on the vertical 
aerosol optical depth obtained with laser shots from the 
location of both detectors. The cuts on the remaining 
events used the average atmospheric conditions and the 
records made during the measurement. The rest of the 
cuts address the reconstruction quality. They include a 
minimum viewing angle of more than 20° in both detec- 
tors, more than 5° difference in the shower detector plane 
(to decrease the effect of scattered Cherenkov light), of 
the global flt of less than 15 p.d.f and bracketing of Xmax 
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within the observed tracks. The application of the same 
cuts to a Monte Carlo data set gives a Xmax resolution 
of 30 g/cm^ and energy resolution of 13%. 

This analysis was presented in a combin ation with an 
earlier result f( Abu-Zavvad et al . 2000bl) ) obtained by 
the HiRes prototype working in coincidence with the 
MIA air shower array. The elongation rate between 10^^ 
eV and 10^^'^ eV was measured to be 93±8.5 g/cm^ with 
a systematic uncertainty of 10.5 g/cm^. This result sug- 
gested a quick transition from heavy to light cosmic ray 
composition. 

The HiRes 2005 paper measured Diq — 54.5±6.5 
g/cm^ consistent with the values for constant composi- 
tion from different hadronic interaction models. The data 
points above 10^^ eV agreed within the errors with the 
results from the HiRes-MIA coincidence experiment. The 
HiRes points, which are derived from fitting the shower 
profile with the Gaisser-Hillas formula, are shown with 
empty circles in Fig. 1191 The lines in the same figure show 
the expectations from three different hadronic interaction 
models: EPOS 1.99, QGSjet II, and SIBYLL 2.1. The 
fraction of protons is different for these interaction mod- 
els. QGSjet II shows almost pure proton composition 
with a possible small He contamination. In the case of 
SIBYLL 2.1 the fraction of protons is smaller but still sig- 
nificant. QGSjet II has a moderate cross section energy 
dependence and fast multiplicity increase. SIBYLL 2.1 
has a fast cross section growth and a relatively low multi- 
plicity. The comparison with EPOS 1.99, which has the 
largest Z^io, may even point at large, but decreasing with 
energy, fraction of protons. 

The interpretation of the X^ax data will become better 
after the LHC results are accounted for in the hadronic 
interaction models used in the analysis. The current 
models do not disagree with each other in the energy 
range studied in accelerators. After the normalization of 
these models to the LHC data the differences in the inter- 
pretation of the experimental results will significantly de- 
crease. An extensive description of the current hadronic 
interaction models used for air shower analysis can be 



found in (jEngel et all 120111 ) 



HiRes also studied the width of the X„iax distributions 
in different energy bins. Proton showers do have a wider 
distribution while iron showers have widths lower by at 
least a factor of two. The comparisons with simulations 
using QGSjet 01 were consistent with a proton fraction 
of 80% and those with SIBYLL 2.1 suggested a proton 
fraction of about 60%. EPOS 1.99 did not exist at that 
time. The general conclusion of the HiRes X^ax study 
is that the cosmic ray composition was heavy at 10^^ 
eV, progressed to light one in one order of magnitude in 
energy and stayed light with a proton fraction from 60 
to 80% above lO^^ eV. 

Although the Auger Collaboration has presented 



its X„ 



measurements at different conferences the 



first journal paper on this topic appeared in 2010 



f (| Abraham et al. . 2010cf) V The results of this study are 
unfortunately very different from those of HiRes. This 
analysis is done using only hybrid events, i.e. events de- 
tected by one or more fiuorescent telescopes plus at least 
one surface detector. Using the timing of the surface 
detector vastly improves the quality of the reconstruc- 
tion. The light collected by the fluorescent detector is 
corrected for attenuation using the atmospheric monitor- 
ing devices. The fitting is done using the Gaisser-Hillas 
function. Events with light emission angle less than 20° 
are not used. Neither are events where Xmax uncertainty 
due to shower geometry and atmospheric conditions is 
more than 40 g/cm^. The limit on the reconstruction 
is set to less than 2.5 p.d.f. The resulting Xmax resolu- 
tion above several EeV is about 20 g/cm^. This number 
is consistent with the checks with stereo fluorescent de- 
tector observations. 

Using data taken between 2004 and March 2009 there 
are 3754 events above 10^* passing all cuts. The highest 
energy event is of energy 59±8xl0^^ eV. The measured 
Xmax in 10 logarithmic bins per decade are shown in 
Fig. [19] with full squares. The elongation rate of the 
three points below lO^^-^s eV is 106^21 g/cm^ and that 
above this point is 24±3 g/cm^. Both these values are 
determined with good fits. Systematic uncertainty is 
around 10 g/ cm^. In absolute value this data set does not 
a ppear extremely diffe rent from the HiRes 2005 analysis 
(^ kbbasi eii"^ . l2005b[ )l but its interpretation is. Instead 
of a constant elongation rate of 54.5 g/cva? we have a 
large one, maybe similar to that of HiRes-MIA, in the 
lower energy part and a short one at higher energy. The 
cosmic ray composition thus have to become lighter up to 
2Q18.25 gY g^j^^ then consistently heavier up to the highest 
energy measured. The Auger Collaboration is careful 
enough to state that such interpretation is reasonable if 
there are no drastic changes in the hadronic interactions 
in this energy interval. 

The Auger collaboration also studied the width of 
the Xmax distributions (RMS) in the same energy bins. 
While the RMS's in the first five bins look consistent with 
a light composition at higher energies there is a steep de- 
crease of T(MB{Xmax) consistent with heavier and heav- 
ier composition as shown in Fig. [20l The RMS values 
decrease from about 55 g/cm^ to 26 g/cm^ in the last 
bin. This distribution width is consistent with cosmic ray 
composition dominated by iron. The interaction model 
predictions for proton showers give 60 g/cm^ with a slight 
energy dependence and these for iron showers give about 
22 g/cm^. The linear decrease of JiMS{Xmax) is not, 
however, consistent with a simple change of the cosmic 
ray composition from pure proton to pure iron. 

It is worth noting that the Auger Collaboration has 
also attempted to use the surf ace detector data for stu d- 
ies similar to those of Xmax ( (| Abraham et al. . 2009dl )). 
Since these results are preliminary we will only give the 
general idea. The width of the shower front depends on 
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the depth of the shower maximum. One can study the 
shower front width by measuring the rise-time of the sur- 
face detector signals. The attempt to do that is fully con- 
sistent with the more detailed fluorescent detector anal- 
ysis. 
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FIG. 19 Depth of maximum measurements of UHECR by the 
HiRes collaboration (2005) analysis shown with empty circles 
and 2010 analysis with full circles and the Auger collabora- 
tion - full squares are compared with the predictions of three 
different interaction models for H and Fe. 

The final analysis of the stereo measure ments of HiRes 



in th e period 1999-2006 was published in (jAbbasi et al 
2010l ). The cuts on the data are more stringent than in 
the previous analysis. Apart from the good weather re- 
quirement, they limit the chance of noise coincidence to 
less than 1% and the longitudinal development fit fo 
less than 4 p.d.f. The final data set of 815 events includes 
only events with zenith angle uncertainty of less than 2°, 
Xmax uncertainty of less than 40 g/cm^, zenith angle 
less than 70° and distance to HiRes II more than 10 km. 
The measured Xmax should be bracketed by the HiRes II 
field of view and have shower detector plane between 40° 
and 130°. The application of the vertical aerosol opti- 
cal depth hourly measurements to the amount of light 
received by the detectors requires a mean upward correc- 
tion of ~15% to shower energy for an event 25 km distant 
from the observatory. Shower segments with emission an- 
gles of less than 5° of a bin pointing direction are not used 
in the analysis. 

The measured light profile of the shower is fit to a 
Gaussian function of the age parameter s to determine 
the shower energy and Xmax- The claim is that the use 
of the Gaisser-Hillas function does not change the results 
within the errors. Showers of energy between 1.6x10^^ 
eV and 6.3x10^^ eV are included in the analysis. 

All uncertainties in the X^ax measurement come from 
the treatment of simulated showers after the detector is 
accounted for. Comparisons of the reconstructed X„iax 
with the original one showed that the selection and recon- 
struction results in Xmax shallower by about 15 g/cm^ 



than the original one. For this reason for the interpreta- 
tion of the measurements the predictions are appropriatly 
scaled. The Monte Carlo measured uncertainty of Xmax 
is better than 25 g/cm^ over most of the energy range. 

This analysis finds a constant elongation rate of 47.9±6 
g/cm^ with fit of 0.86 p.d.f over the whole range with 
systematic uncertainty of 3.2 g/cm^. Most of the system- 
atic uncertainty is due to the event selection cuts. 

HiRes also presents the energy dependence of the Xmax 
fluctuations in the same energy bins. These numbers 
are obtained in a different way from those of the Auger 
Collaboration. Since these Xmax distributions are wide 
and asymmetric, the HiRes analysis fits them to Gaussian 
distributions truncated at 2xRMS. The distributions are 
still wide as shown in Fig. [20l 
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FIG. 20 Width of the Xmax distributions as measured by 
HiRes (full circles) and Auger (full squares). Note that the 
width is presented in different ways (see text) and the points 
cannot be compared directly. 



The heavy cosmic ray composition derived from the 
Auger data suggests that the strong decline of the cos- 
mic ray flux may be caused by exceeding the maximum 
acceleration energy at the cosmic ray sources. In such 
a case only iron nuclei could be accelerated to energies 
exceeding 10^° eV. 



D. Transition from galactic to extragalactic cosmic rays 

One of the reasons for identifying different features at 
the end of the cosmic ray spectrum is to study the transi- 
tion between the galactic and extragalactic components. 
The common opinion is that most of the cosmic rays 
above 10^^ eV are of extragalactic origin and the GZK 
feature supports that. The main question was (and is) 
the origin of the dip at around 3x 10^^ eV. The prevailing 
school of thought was that the dip is at the intersection of 
the Kala.ctic a nd ext ragalactic components as exp lained in 
(|Hillasl . [l983) and (iBahcaU and Waxmanl . l2003l ). In this 
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model the extragalactic cosmic rays have a flat spec- 
trum and the galactic ones have a steep E~^-^ spectrum 
as shown in the upper panel of Fig. with two values 
(3&4) of the parameter m in the evolution described as 
(l + z)™ up to redshift of about 2. As is seen from the fig- 
ure the cosmological evolution does not affect much the 
predicted spectra because the observed UHECR have to 
be local. The galactic cosmic rays, although with a small 
contribution, extend well above 10^^ eV as shown with a 
dashed line in the figure. It became soon obvious that ex- 
tragalactic cosmic rays cannot have such a flat injection 
spectrum subsequent models deal with E"^ '^ or slightly 
steeper spectra. 



10 = 



10 = 



(0 
(0 

CM 

i. 1023 
UJ 

z 

■o 
m 
UJ 




10^ 



10 



y=1, m=3,4 



^o8oo 



I I I I 



y=1.7, m=0 




23 

10 



J I 



J I 



18 



10 



19 



1020 



10 



21 



E(eV) 



FIG. 21 The mod e ls of jBahcall and WaxmarJ, I2003D (a) and 
(jBerezinskv et all l2005f ) (b) compared to the more recent 
data of Auger (full squares) and HiRcs ( circles). 



Soon after that (jBerezinskv et all 120051 ) suggested 



a totally different model. The dip is caused by the 
pair production interaction of the extragalactic protons 
with CMB as predicted by (jBerezinskv and Grigor'eva , 
19881 ). The model u nderwent some development 
later fl Aloisio et al . 200?! )) to convey some of its details. 
The shocking part of this model, illustrated in the lower 
panel of Fig. [21] is that the cosmic ray injection spec- 
trum was a steep E~^-^ rather than the expected flat 
one. At about 10^^ eV the injection spectrum has to be- 
come much flatter for the flux not to exceed the measured 
cosmic ray spectrum. The transition from galactic to ex- 



tragalactic cosmic rays should then happen below 10^^ 
eV. There is no need for cosmological evolution of the 
cosmic ray sources in the model. Extragalactic cosmic 
rays had to be almost exclusively protons. 

At about the same time a third model for the tran- 



sition became available ( (jAUard et all [20071)) following 
the calculations of the heavy nuclei propagation in extra- 
galactic space and th eir interactions wit h the CMB and 



other photon flelds f ijAUard et aLl . l2005l ): (jHooper et al. 
20071 )). Since the extragalactic cosmic rays in such a 



model may have at least five chemical components, i.e. 
many more parameters, the model is much more compli- 
cated, although it could be made to fit the cosmic ray 
spectra as well as the first two. These propagation cal- 
culations also showed that protons and iron nuclei have 
approximately equal energy loss lengths, while all inter- 
mediate nuclei would disintegrate at much shorter dis- 
tances. 

Since the spectrum shape alone cannot answer the 
questions about the transition, the answer could only 
come from accompanying composition study. The chem- 
ical composition derived by the Hi Res experiment to- 
gether with the measurement of (lAbu-Zavvad et al. 



2000bl ) would claim that the galactic cosmic ray spec- 
trum does not extend above 10^^ eV and higher energies 
contain only protons with a small admixture of light nu- 
clei. This admixture may be different in the 2005 and 
2010 analyses but it seems to be constant and belong to 
the same population. Such interpretation may not be 
consistent with a proton cosmic rays knee at 3x10^^ eV 
or lower as derived by the Kascade experiment since in 
such a case the iron knee would start at 10^^ eV and 
would leave about one o rder of magnitude of energy un- 
explained ( (jHillasl . [2005l )). 



A similar simple interpretation of the Auger Xmax and 
RMS result is impossible. The elongation rate derived 
from the first three points seems to show a quick tran- 
sition from heavy to light nuclei followed by a slower 
transition to heavy nuclear composition. 

We first have to understand if at least a part of the 
Xmax behavior is not due to a sudden change of hadronic 
interactions at close to 50 TeV, well above the LHC 
maximum energy. Then we have to relate the change 
of the injection composition to the shape of the energy 
spectrum that is in this range quite different from that 
of HiRes. These are not simple problems to solve and 
in our opinion they will take years. The main hope is 
that Auger and HiRes would examine each others X^ax 
analysis techniques and will come up with similar, if not 
identical, X^ax values as a function of the energy. The 
lower energy extensions of the UHECR arrays described 
in Section 1x1 and the use of different composition related 
parameters (muon to electron density ratios) may also be 
of big help. 
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VIII. SEARCH FOR THE SOURCES OF UHECR 

Before describing the searches for the sources of 
UHECR we will briefly introduce some of the ideas about 
the strengths of the magnetic fields in the Universe. 
These are important parameters because particle scat- 
tering in the magnetic fields can hide the sources. A 
nice review of the investigations and results of the stud- 
ies of astrophysical magnetic fields can be found in (jBeckl , 
200 ih . 



netic field above and below the galactic plane. More re- 
cent measurements tend to show an extended halo that 
c an contribute a lo t to the cosmic ray scattering angle 
( (I Jiang et aLl . |2010[) ). A standard way to study UHECR 



A. Galactic magnetic fields 



scattering angle is to inject negatively charged nuclei in 
a magnetic field model a nd follow their trajectories until 
they leave the Galaxy fi Stanevl . Il997t )). Such exercises 
with different toy galactic field models give scattering an- 
gles at 100 EeV between 2° and 4° depending on the cos- 
mic ray direction. Some other estimates, however, give 
much higher values, up to 10° {R. Beck, private commu- 
nication). 



Galactic magnetic fields are very important for the 
scattering of UHECR because they definitely have a large 
scale structure. This means that cosmic rays coming from 
the same direction will scatter in a similar way and the 
scattering will shift the arrival direction away from the 
true source. 

The regular magnetic field strength is measured mostly 
by studies of the Faraday rotation of the radio emission of 
pulsars. The rotation measure RM , measured in rad/m^ 
is proportional to UeB^^dl where d is the distance to 
the source and n,, is the electron density. The integral 
over the measured or assumed Ue is used to extract the 
magnetic field strength as de scribed in a recen t review of 
the galactic magnetic field f( Han et al. . 2006[) V In prin- 
ciple the magnetic field strength is proportional to the 
matter density in the Galaxy and is decreasing with the 
galactocentric distance. The field decrease in the galac- 
tic plane is best described with an exponential function 
g-flGc/8-5 T^Y^QTf-Q |;jjg distance from the galactic center 
Rgc is in kpc. The local regular field in the vicinity 
of the Solar system has a strength of about 2 /xGauss 
and points counterclockwise close to the direction of the 
Carina-Sagittarius arm. 

This expression is valid for Rgc bigger than 3 kpc be- 
cause the field inside that circle is difficult to study and 
is not well known. In regions near the galactic center 
mGauss fields have been observed pointing almost per- 
pendicular to the galactic plane. This led to suggestions 
that there is a strong magnetic dipole in the galactic cen- 
ter. At the solar system the dipole field strength is about 
0.3 ^Gauss and points North in galactic coordinates. 

The more general estimates of the total field strength 
at our location give values of 6/xGauss which results in 
random field strength twice as big as the regular field. 
There are also ideas that the random field reaches max- 
imum inside the galactic arms (because of the stellar 
fields pointing in different directions) and the regular field 
reaches maximum in the inter-arm space. The random 
field is not very important for UHECR scattering because 
its scale size is only 50-100 pc. 

A very important question which is far from solved is 
the galactic magnetic halo, i.e. the extension of the mag- 



B. Extragalactic magnetic fields 

Our knowledge of the extragalactic magnetic fields is 
much sma ller and still on t he basic level of the excellent 
review of ( KronberS , 1994 ). Although /xGauss magnetic 
fields have been observed in clusters of galaxies, such 
objects enclose a small fraction of the Universe (10~^ or 
less) and the upper limit of the average magnetic fields 
is 10 ~^ Gauss = 1 nGauss if the correlation length of the 
field Lc is 1 Mpc, the average distance between galaxies. 
But even such small fields can affect the propagation of 
UHECR. 

The angular deflection due to random walk 6 would 
then be 

-1/2 



2.5° E-^ B_q d 
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(23) 



where £^20 is the energy in units of 10^° eV, B-q is the 
magnetic field strength in nGauss, dioo is the source 
distance in units of 100 Mpc and Lc is the correlation 
length in Mpc. The random walk causes a propagation 
path length Ad that is larger than the distance to the 
source and causes increased energy loss. It depends on 
the square of the parameters above and is 



Ad 



0.047 i;2o^B? 



9 '^loo Lc Mpc. 



(24) 



The increased propagation distance causes a correspond- 
ing time delay. In case UHECR are generated in a GRB, 
or in an active state of an AGN we may not be able to 
correlate these events with the resulting cosmic rays. 

The situation changes drastically if the UHECR en- 
counters an extended region with organized magnetic 
field. In principle this should be a rare occasion except 
close to a powerful astrophysical system where such fields 
have been observed. Depending on the field strength, its 
direction toward us, and structure of the field, the angu- 
lar deflection could be much larger. 



C. Correlation of the arrival directions of UHECR with 
astrophysical objects. 

The first attempt to correlate the arrival direction of 
UHECR with known astrophysical objects was published 
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by (jStanev et aUll995[ ). The authors used 143 events of 
energy more than 2x10^^ eV detected by the Haverah 
Park array, together with the statistics of the Vulcano 
Ranch, Yakutsk and the prehminary data of AG AS A. 
The authors studied the angular distance between the 
UHECR events and the super-galactic plane (SGP), 
which is the plane of weight of almos t all extragalac- 
tic ob jects within redshifts below 0.04 (( de Vaucouleursl . 
1956I )). The conclusion was that at energy above 4x 10^^ 
eV the average and RMS distances of UHECR to SGP 
are much closer than would be expected from an isotropic 
distribution of the UHECR sources. 

With the increase of the AGASA statistics that started 
to dominate in the late 1990s the correlation with the 
SGP decreased. Other effects were claimed by that 
experiment: a large scale isotropy and small scale 
anisotropy. The anisotropy was defined by the fact 
that three pairs and a triple of events coming withi n 
2.5° of each other were found ( (|Takeda et all Il999l )) 
among 47 events of energy above 4x 10"'^^ eV. The chance 
probability of this happening from isotropic distribution 
was less than 1 %. Soon after this clustering analysis 
was extended to inclu de the previously detected events 
( (jUchihori et allhOOdi} ) . The conclusions from that anal- 
ysis were slightly different. Since the angular resolution 
of the older experiments was worse the clustering was 
analyzed in angular distances of 3, 4, and 5°. Twelve 
doubles and 2 triples were found within 3°. The empha- 
sis, though, was on the fact that 8 of the doubles and the 
2 triples lie within 10° of the super-galactic plane, which 
had a chance probability of 0.1 and 0.2 % for an isotropic 
distribution of the sources. 



1. Correlation of the Auger events with AGN 

After collecting and exposure of 4,390 km^ sr yr the 
Auger collaboration noticed that many of their higher en- 
ergy ev ents are close to the activ e galac tic nuclei from the 
VCV f (|Veron-Cettv and VeronI , |2006|) ) catalog. They 
did a scanning analysis varying the angular distance, 
event energy and the source distances. The best correla- 
tion appeared for angular distance of 3.1°, event energy 
above 5.7x10^^ eV (57 EeV) and distance of 75 Mpc 
(redshift less than 0.018). The same analysis was re- 
peated when the exposure more than doubled and the 
total number of high energy events reached 27. Twenty 
of these events were within 3.1° of the AGN from the 
VCV catalog when only 7.4 were expected for isotropic 
sources. T he chance probability for this h appening was 



1.7x10^3 C ([Abraham et all l2007bl. l2008al )). A signifi- 



cant number of events were close to the nearby (distance 
of 3.8 Mpc) radio galaxy Cen A. There were no events 
close to the powerful AGN M87. When the events with 
galactic latitude less than ±12°, where the catalog cov- 
erage is smaller and UHECR scattering in the galactic 



magnetic field is supposed to be stronger, were excluded 
the strength of the correlation increased and 19 out of 
21 events correlated with at least one AGN as shown in 
Fig. [2H This strong correlation was surprising because 
of several reasons. First of all the VCV optical catalog 
includes many low power objects that are not likely to 
accelerate particles to such high energy. Secondly, the 
0.018 redshift does not correspond to the GZK horizon 
(the distance up to which cosmic ray sources contribute 
significantly to the flux observed above a certain energy) 
for energy of 57 EeV and the question arose if the Auger 
energy scale was low by about 25% which would bring 
the two distances close to each other. 

The Auger collaboration did not claim that the AGN 
from the VCV catalog are the actual sources, which may 
have a sky distribution similar to that of the correlating 
AGN. The anisotropy of the UHECR sources was em- 
phasized in the papers. 

The analysis was r epeated by the HiRes experiment 
fl Abbasi et al . 2008cl )) as close to the original as possi- 
ble. There were only two out of 13 events with similar 
energy that correlated with the same AGNs and the con- 
clusion was the opposite. The HiRes field of view is not 
the same as that of Auger and the VCV catalog has dif- 
ferent coverage of the corresponding fields of view. Still 
the results from the two analyses appeared to be con- 
troversial since HiRes sees one half of the Auger field of 
view. 

Recently the Auger collaboration presented the cor- 
r elation results from an exposure of 20,370 km^.sr.yr 
C dAbreu et ali I2OIOI) ') which contains 69 events of en- 
ergy above 55 EeV (corresponding in the contemporary 
energy assignment to 57 EeV in 2007). The complete 
catalog of the 69 Auger events published to date are 
shown in Fig. [511 The correlation is now weaker - 42% 
of all events (29/69) correlate, compared with 74% in 
2007. The event reconstruction is constantly improving 
and because of that a couple of events move from one 
group (above 55 EeV or not) to the other. If the events 
participating in the initial parameter scan are excluded 
the corresponding fractions of correlated events are 69% 
and 38% respectively. This does not mean, though, that 
the observed UHECR are isotropically distributed as the 
expected fraction of correlated events is 21%. 

It should be noted that the scattering in the extra- 
galactic magnetic fields sho uld be much stron ger accord- 
ing to some calculations. ( Rvu et al. . 2010l ) predict an 
average scattering angle of 15° for cosmic ray protons 
above 60 EeV. 



2. Correlation with sources from other catalogs 

The most difficult part of the search for UHECR 
sources is that we have no idea what is the best proxy 
for cosmic ray acceleration to the highest energy - is 
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FIG. 22 Correlation of the arrival directions of UHECR with AGN from the VCV catalog. The shaded part of the sky is not 
visible by Auger. The gray squares are the AGN within z less than 0.018. The Auger events are shown with circles. The first 
27 events are half filled . The 13 HiRes events are shown with black dots. The thin lines show the six regions of the sky to 
which Auger has equal exposure. The wide gray line is the supergalactic plane. 



it the optical/UV luminosity, or the X-ray one, or still 
higher energy 7-ray emissivity? The last paper on the 
Auger events anisotropy explores the cor relations with 
two more catalogs: the 2MRS catalog f ijHuchra et al\ . 
20051 )). which contains the brightest galaxies from the 



2 MASS cata log, and the Pal e rmo S wift-BAT hard X- 
ray catalog (( Cusumano et aA . I2OIOI )). 2MRS contains 
13,000 galaxies within 100 Mpc and 22,000 galaxies 
within 200 Mpc. The Swift-BAT catalog has the advan- 
tage to cover well the region of the galactic plane. It con- 
tains 133 extragalactic sources within 100 Mpc and 267 
within 200 Mpc. The correlation of the Auger UHECR 
arrival directions with the positions of the objects in both 
catalogs is much better than an isotropic source distri- 
bution would suggest. 

To fully understand the correlations with catalogs con- 
taining different number of objects and to estimate the 
statistical significance of these correlations the Auger col- 
laboration used a different approach. The catalogs were 
used to create maps of possible sources using the object 
densities per unit area of sky where each object position 
was extended by several degrees. These extensions are 
supposed to account for the particle scattering in mag- 
netic fields and the angular sensitivity of the experiment. 
The UHECR luminosity of the sources were scaled with 
the distance and with the observed source luminosity at 
different wavelengths. With the use of simulations the 
events were then separated in source and isotropic frac- 
tions with different confidence levels. The isotropic frac- 
tion became on the average 0.64 for the 2MRS catalog 
and 0.62 for Swift-BAT with huge error bars even at la 
level. In a way this analysis produced similar results to 
the contemporary correlation with the VCV catalog. 



The last test of isotropy was made with studies of self 
correlation - a comparison of the number of event pairs 
as a function between the angular distance of the pair 
compared to that of isotropic source distribution. The 
number of experimental pair events is consistently above 
the expectations. The biggest deviation is at angular 
distance of 11°, where the experimental events show 51 
pairs while 34.8 pairs are expected for an isotropic distri- 
bution. At angular distances higher than 45° the number 
of pairs is consistent with isotropy but below 30° it is not. 



3. Events coming from specific objects 

Ever since the publication of the first analysis of the 
correlation of the Auger events with extragalactic objects 
the question was why there are so many events coming 
from directions close to Cen A and there are none from 
the Virgo cluster and M87. With the increased contem- 
porary statistics Auger was able to analyze better this 
fact. There are still no events coming from less than 18° 
from M87. And there are now 13 events coming from 
less than 18° from Cen A and two events very close to it. 
Cen A is close to the direction of the Centaurus cluster 
but is not a part of it. 

M87 is almost 5 times more distant than Cen A, which 
is at a distance of 3.8 Mpc. It is also in a region where the 
Auger exposure is three times less as shown in Fig. [22] 
Using these two rough numbers one expects 75 times less 
events from M87 than from Cen A. In other words one 
expects 13/75 events coming from M87 if it has the same 
CR luminosity as Cen A. The lack of events then is not 
a problem. 
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upper bound of the ultrahigh energy muon neutrinos 
and antineutrinos would be if UHECR lose all of their 
energy to neutrino production. Using the average energy 
loss they arrived at an upper bound of 



EldNy/dE^ = 1.5 X 10"^GeV(cm2.s.sr)" 



(25) 



which after accounting for the cosmological evolution of 
the sources as z)^ is increased by a factor of 3. The 



limit was criticized by (jMannheim et all l200lf ) who de- 



rived a more realistic limit that only touched the limit 
of Waxman&Bahcall at 10^^ eV. Both limits are shown 
with thick gray lines in Fig. [24l 



FIG. 23 Number events as a function of the angular distance 
to Cen A. The thick gray line shows the expectation from 
isotropic distribution of the cosmic ray sources. 



The 13 events coming from directions close to Cen A 
are mostly responsible for the excess of self correlation 
discussed above. The events coming from this direction 
have 28 pairs coming with separation less than 11°. For 
an isotropic distribution one expects 3.2 events rather 
than 13, while the map based on the 2MRS catalog pre- 
dicts 9.2 and that based on Swift-BAT catalog predicts 
20.6. 

Figure [23] shows the comparison of the number of 
events coming at different distances from this object com- 
pared to the expectations from an isotropic distribution. 
A Kolmogorov-Smirnov test of this distributions estab- 
lishes 96% significance or about 2a deviation from an 
isotropic distribution. The question then is if at least 
a part of these events come from the Centaurus cluster 
rather than from Cen A. This does not appear likely be- 
cause the Centaurus cluster is further away than Virgo 
and one would expect a small fraction of events coming 
from there for equal CR luminosities, which is of course 
not guaranteed. 



IX. ULTRAHIGH ENERGY NEUTRINOS 



The 
trahigh 

by (jBeresinskv and Zatsepin . 



relationship between 
energy neutri nos was 



important relation between 
flux a nd the flux of di ff use 
rived (IWaxman and Bahcalll 1999() 



UHECR and ul- 
first noted 
19691) . Later on, an 
the observed UHECR 
neutrino was de- 
The authors 



( Waxman and Bahcall . 
took a simple basic approach to the problem and using 
the measured cosmic ray flux at 10^^ eV. Then they 
assumed a flat, 7=2 injection spectrum and calculated 
the emissivity of UHECR in the Universe which came to 
10^* erg/(Mpc3.yr) in the range lO^^ - lO^i eV. The next 
observation is that a fraction of these cosmic rays would 
have photo-production interaction at their sources and a 
fraction of their energy loss e would go to neutrinos. The 



A. Cosmogenic neutrinos 

Cosmogenic neutrinos wer e flrst suggested 
by (jBeresinskv and Zatsepiril . Il969h . These are neutrinos 
that are produced by UHECR in photo-production 
interactions in the CMB and other photon fields in 
propagation. The original paper did not produce much 
interest since the contemporary experiments could 
not detect high energy neutrinos. The shapes of the 
cosmogenic neutrino spectra are very different from 
those of the Waxman&Bahcall limit. Muon neutrino 
and antineutrino spectra peak at about 10^^ cV and 
significantly decline at both lower and higher energy. 
These spectra are shown in Fig [24] together with the 
two limits assuming the same astrophysical input. Even 
after the multiplication by E^ the electron neutrino and 
antineutrino spectra show an extension to lower energy 
which is due to Df, from neutron decay. 
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FIG. 24 The upper lim its on the ultrahigh e nergy neu- 
trino fluxes derived by (iWaxman and Bahcalll . Il999f ) (la- 
belled) and ([Mannheim et al\ . l200ll ) are compared to a cal- 
cul a-tion of cosmogeni c neutrinos produced by UHE protons 
by ([Engel et aZ.I. l200ll ). Electron neutrinos and antineutrino 
fluxes are plotted with a dashed line. One can see the contri- 
bution of neutron decay at lower energy. 

The magnitude of the cosmogenic neutrino spectra de- 
pends on the cosmic ray injection spectrum and compo- 
sition, on the distribution of UHECR sources, and very 
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strongly on the cosmological evolution of these sources. 
For flatter injection spectra more UHECR can undergo 
photo-production interactions and hence generate more 
neutrinos. Cosmological evolution importance has a sim- 
ple explanation - UHECR can arrive to us only from very 
low rcdshifts (less than 0.05) while neutrinos can travel 
without energy loss (except adiabatic) from the whole 
Universe. If the cosmological evolution of the UHECR 
sources peaks at z=2, as it does in many models, the 
cosmogenic neutrino production would peak close to z=3, 
when the source emission is much stronger. 

The influence of the cosmic ray composition on the 
cosmogenic neutrino flux is even stronger, although more 
difficult to evaluate. Fig. shows the fluxes of cosmo- 
genic neutrinos calculated for UHE protons. The solid 
line shows the sum of muon neutrinos and antineutrinos 
and electron neutrinos and the dash dot line shows the 
flux of electron antineutrinos from neutron decay. The 
input parameters come from the Auger energy spectrum 
flt that produces the larger amount of cosmogenic neu- 
trinos (protons, 7=1.3, cosmological evolution (1 + z)^). 
If UHECR are not all protons, the solid line should come 
down keeping all other parameters stable. For 20% pro- 
tons in UHECR the flux would be lower by a factor of 
5. At the same time the flux of cosmogenic would rise 
since heavy nuclei photo-disintegrate and emit neutrons 
that decay. The estimate of the increase of the Pe flux is 
more complicated but it would increase roughly also by 
a factor of 5. The cosmogenic neutrino flux would then 
be dominated at production by this neutrino flavor. Af- 
ter propagation it would be shared equally by all three 
neutrino flavors. 

10"^ 1^ -^r- ^T-^ --T-^ --T-^ --T^ -^r- ^ 
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FIG. 25 Neutrino fluxes calculated from the Auger energy 
spectrum interpretation that produces the larger amount of 
cosmogenic neutrinos. The solid line shows the sum of i^,, 
+ Vfi + Ve- The dash-dot line shows the flux of i/e- Any 
heavier nuclei contribution to the UHECR flux can only raise 
the dash-dot line and decrease the solid one. The thick gray 
line shows the median (halfway between optimistic and pes- 
simistic) limit set by Auger on Ut- 

The conclusion is that a measurement of the cosmo- 



genic neutrino flux is a complimentary measurement to 
that of the UHECR spectrum and composition. Even the 
detection of a few events will considerably help the analy- 
sis of the UHECR features and origin. The problem here 
is that the current limits on the UHE neutrino fluxes are 
generally above the predictions shown in Fig. [55] The 
measurement of the cosmogenic De spectra and their os- 
cillations is even more difficult because of the energy de- 
pendence of the neutrino-nucleon cross section. 



X. REMAINING PROBLEMS AND EXPECTATIONS 
FROM FUTURE EXPERIMENTS 

A. Remaining problems 

It is obvious from the controversial results on the chem- 
ical composition of UHECR that this is the main un- 
solved problem. The uncertainty of the chemical com- 
position also affects the interpretation of the end of the 
cosmic ray energy spectrum. If the UHECR composi- 
tion is dominated by protons the most likely explanation 
is the GZK effect. If, however, the composition is in- 
creasingly heavier to the highest energies it could be a 
result of reaching the maximum acceleration rigidity at 
the UHECR sources. 

At the highest energies (above a few tens of EeV) two 
important observables need to be measured with better 
precision, the composition and the anisotropics. Both 
will tell us about the sources and their distribution as well 
as about the mechanisms at play in accelerating particles 
up to 100 EeV or more. 

Currently the available data from Auger and HiRes 
appear contradictory and no model is able to explain in 
a coherent way all the observations. Moreover, in the 
Auger data Centaurus A is today the sole possible source 
candidate that may have been seen in the sky. Can this 
possibility be confirmed, is CenA the only source visible 
from the Southern hemisphere? From both hemispheres? 

At the highest energies the effort needs to be pursued 
along at least three lines: covering the whole sky, in- 
creasing the statistics by instrumenting larger surfaces 
or volumes, and improving the measurements adding 
new detector components. To make definite progress, 
the next generation of detectors should be able to mea- 
sure independently, and if possible redundantly, all EAS 
components. This includes in particular, electromagnetic 
shower profile with a maximum of a few tens of g/cm^ res- 
olution, as well as the muonic and electromagnetic com- 
ponents at ground to better constrain hadronic model 
and the first interaction dynamics. 

At the EeV scale, the expected transition from galac- 
tic to extra-galactic origin in the cosmic ray spectrum 
has not been confirmed. Several features in the energy 
spectrum need attention. Is there a second knee around 
0.1 EeV ? or at almost 1 EeV as measured by the Akeno 
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array (jNagano et all Il992f) ? How pronounced is the 
ankle ? What is its origin? Today the interpretations 
in terms of a pure proton composition undergoing e+e~ 
pair creations, or in terms of the galactic to extra galactic 
transition of a mixed composition seem equally (in)valid. 
What is the level of anisotropy in this energy range, can 
the above models a ccommodate or predict the low values 
already reported (( Abraham et al . 2009bl )). Again the 



only hope for light can come from more accurate mea- 
surements of this regions, both in terms of statistics and 
in terms of multi-parametric observations. 



B. Extensions of Auger South 

The High Eleva tion Auger Telescope (HEAT, 
(|Kleifges et all l2009t )) and the A uger Muon and Infi ll 
for the Ground Array (AMIGA, (jPlatino et all l2009t )) 
have been added to the original design of the Pierre 
Auger Observatory. Improving the efficiency of the ob- 
servatory in the 0.1 to 1 EeV range, these extensions will 
efficiently test the various models for the acceleration 
and transport of galactic and extragalactic cosmic rays 
in the transition region. 

Studies of this region require not only a better collec- 
tion efficiency to improve the statistics but also powerful 
mass discrimination capabilities. While very high energy 
showers can be efficiently measured by fluorescent tele- 
scopes from distances up to several tens of kilometers, 
lower energy ones do not emit sufficient light to bee seen 
further than a few kilometers away. For the same position 
of Xmax closer showers appear at higher elevation than 
distant ones. Since low energy showers reach their maxi- 
mum of development faster, coverage above the 30° limit 
of the original Auger telescopes is required. HEAT is 
composed of three fluorescence telescopes of the same ba- 
sic design as the original Auger telescope and is installed 
at the western fluorescence detector site (Coihueco) of 
the observatory. They can operate in two positions. Hor- 
izontally they share the same field of view as the original 
telescopes. This position is used for laser and drum cali- 
bration of the instruments as well as for inter-calibration 
using shower data. Tilted upward by 29°, this is the 
normal operation mode in which the nearby upper part 
of the atmosphere is observed. Construction took place 
in 2008-2009 and first light was seen from one of those 
telescopes in January 2009. 

Routine observation with the HEAT telescope began 
in 2010. Figure [26] shows the longitudinal shower profile 
of an event recorded in coincidence with the Coihueco 
telescope. The reconstruction of this event gives a shower 
energy of 0.2±0.02 EeV and a distance of 2.83±0.06 km 
from Coihueco. It is clear from the plot that the data 
points provided by the HEAT telescope are mandatory 
to properly reconstruct the shower development profile. 

The Auger observatory reconstruction is based on the 
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FIG. 26 Longitudinal profile of a 0.2±0.02 EeV shower re- 
coded in coincidence by the HEAT and Coihueco fluorescence 
telescope of the Auger observatory. 



hybrid technique. To provide the HEAT telescope with 
adequate information from the surface array it was nec- 
essary to also increase the surface detector density at the 
foot of the telescope. An infill array of 85 detectors is 
deployed on two grids of one half (750 m) and one fourth 
(433 m) of the regular Auger surface array grid. Mea- 
suring the muon densities on the ground together with 
the electromagnetic component provides important in- 
formation on the cosmic ray composition in addition to 
the longitudinal shower development. Such a multi para- 
metric measurement allows to study independently the 
evolution of Xmax and of the muon densities which are 
linked in a similar way in all interaction models. The 
AMIGA extension aims to provide such information by 
measuring the shower muons with buried muon counters. 
Each counter is made of a segmented plastic scintillator 
read out by wave shifting fibers connected to a 64 chan- 
nels multi anode PMT. 

The muon lateral distribution function is adjusted to 
the counter data to provide the number of muons at 600 
m from the shower axis. Realistic Monte Carlo analysis 
together with an improved reconstruction showed a rela- 
tive precision on the estimated muon density better than 
20% in the energy ran ge of 0.4 to 3 EeV acc essible to 
the 750 m infiU alone f (|Supanitskv et aLl . lioosl ) 1 . 

At the time of writing nearly all of the 750 m infill grid 
is completed and is operating while a muon counter has 
been buried and successfully tested. Completion of this 
effort is expected to take place in 2011-2012. 

Finally, co-located with the infill array, the Auger 
collaboration is currently installing the first phase 
of the Auger Engineer ing Radio Array (AERA, 



van den Berg et al . 20091 )). The base line parameters 



for AERA comprise about 150 radio detection stations 
distributed over an area of 20 km^. The main scien- 
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tific goals of the project are the thorough investigation 
of the radio emission from an air shower at the highest 
energies, the exploration of the capability of the radio- 
detection technique, and the provision of additional ob- 
servables (calorimetric energy and shower profile deter- 
mination with 100% duty cycle) for the composition mea- 
surements between 10^^'^ and 10^^'^ eV 

In order to increase the amount of data on the shower 
longitudinal development, today severely limited by the 
10% duty cycle of the fluorescence detector, the Auger 
collaboration is pursuing several R&D programs aiming 
at measuring the shower longit udinal development us- 



ing microwave radio techniques ('(IGambetta et al. 
Gorham et all . l2008bt IPrivitera et g/J . l2010[ l). 
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C. Telescope Array 

The Telescope Array (TA) is a new hybrid detector 
that started collecting data in 2009 in Utah, USA, at 
39°N, 120°W and ahitude of 1500 m a.s.l. Its surface ar- 
ray (SD) currently consists of 607 scintillator counters on 
a square grid with dimension of 1.2 km. Each scintillator 
detector consists of two layers of thickness 1.2 cm and 
area of 3 m^ . The phototube of each layer is connected 
to the scintillator via 96 wavelength shifting fibers which 
make the response of the scintillator more uniform. Each 
such station is powered by a solar panel that charges a 
lead-acid battery. The total area of the surface array 
is 762 km^. The surface array is divided in three parts 
that communicate with three control towers where the 
waveforms are digitized and triggers are produced. Each 
second the tower collects the recorded signals from all 
stations and a trigger is produced when three adjacent 
stations coincide within 8 /Ltsec.The SD reaches a full ef- 
ficiency at 10^^-^ eV for sho wers with zenith angle less 
than 45 ° f (|Nonaka et al\ .\m)9)). This angle corresponds 
to SD acceptance of 1,600 km^sr. 
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The fluorescence detector (FD) consists of three fluo- 
rescence stations as shown in Fig. [27l Two of them are 
new and consists of 12 telescopes with field of view from 
elevations of 3° to 31°. The total horizontal field of view 
of each station is 108°. Each telescope has a camera 
consisting of 256 PMT with field of view l°xl°. The 
signals are digitized by 40 MHZ sample FADC and the 
waveforms are recorded when signals are found in 5 adja- 
cent PMTs. The third station has 14 telescopes that use 
cameras and electronics from HiRes-I and mirrors from 
HiRes-II. The fluorescent telescopes are calibrated with 
N2 lasers, Xe flashers , and an electron linear accelerator 
r (|Tokuno et alihOO^ ). 

The atmosphere is monitored for clouds by IR cameras 
and with the use of the central laser facility which is in the 
center of the array at 20.85 km from each station. The 
fluorescent stations are positioned in such a way that they 
cover the whole area of the surface detector. The mono 
acceptance of the FD is 1,830 km^sr and the stereo one 
is 1040 km^sr. The total energy resolution is 25% and 
the Xfnax resolution is 17 g/cm^. 

TALE is the lower energy extension of the Telescope 
Array which will be deployed to study cosmic rays of 
energy 10^^'^ to 10^^ eV. It consists of an infill array and 
a fourth fluorescent station inside TA. The field of view 
of this station will be elevations of 33° to 71° so it will 
be able to see X^ax of the lower energy showers. The 
infill will consist of 100 scintillator counters on a square 
grid of 400 m and muon counters. See its description at 
http:/ /wow. telescope array.org. 



FIG. 27 Sketch of the Telescope array geometry. The surface 
detectors are indicated with full diamonds and the telescope 
stations with arcs. 



D. Auger North 

Based on the same detection principle as the south- 
ern observatory the design of the northern site of Auger, 
or Auger North for short, i s focussed on collecti ng sig- 
nificantly larger statistics (( Bluemer et al . 2010l) ). Its 
target energy lies above the 60 EeV threshold where 
anisotropy in the distribution of cosmic ray sources 
within the GZK sphere has been detected. Motivation 
for such a detector are plenty, they principally concern 
the determination of the cosmic ray composition up to at 
least 100 EeV. This would lead, on one hand, to the iden- 
tification of the trans-GZK cosmic ray sources and their 
acceleration mechanisms and, on the other hand, to the 
study of particle interactions at center of mass energies 
far beyond any man made accelerators. 

Covering 20,000 km^ this new facility is to be de- 
ployed in the state of Colorado (USA) in the north- 
ern hemisphere to provide the Auger collaboration with 
full sky coverage. Composed of a particle array of 4000 
Cherenkov tanks principally on a miles grid and 39 
fluorescence telescope overlooking the atmosphere over 
the whole array. This configuration should reach 90% 
efficiency above 30 EeV for proton primaries. 
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The expected performances of this detector are similar 
to its southern counterpart but at a higher energy due to 
the larger spacing. For example the angular resolution 
is expected to be better than 2.2° above 50 EeV. The 
statistics above 60 EeV is expected to be of the order of 
150 events per year. Out of those, of order of 10 per year 
should have an appropriate profile reconstruction from 
the FD telescopes for mass composition measurements. 

Construction of this facility was planned for 2011 and 
should last 5 years. However as of the end of 2010, the 
funding situation and prioritization in the USA does not 
allow for such construction to start in the short term. 
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E. EUSO - JEM-EUSO 

The Extreme Universe Space Observatory (EUSO) is 
an UV telescope mounted on an external facility of the In- 
ternational Space Station (ISS) which observes the atmo- 
sphere to detect light signals from UHE cosmic rays and 
neutrinos. It is a monocular telescope that measures the 
air shower fluorescence light and the Cherenkov light dif- 
fusively reflected from the surface of the Earth. The ini- 
tial idea of such experiment was proposed and developed 
by th e pioneer of the UHECR research John (Jjinsley, 
19981 ). He was later joined by Livio Scarsi and a group of 
scientists from the University of Palermo. Initially EUSO 
was approved by the European Space Agenc y for a Phase 
A conceptual study (( Catalano et al\ . \200(h ). It was not 
approved to be mounted on the European research mod- 
ule of ISS and was taken over by the Japanese Space 
Agency. It is now known as JEM-EUSO. 

EUSO is a wide angle (±30°) camera with a diameter 
of 2.5 m. The UV light is imaged through Fresnel lens 
optics and detected by a segmented focal surface detec- 
tor using multi-anod PMT. The aim is to have 1 km^ 
resolution on the surface of the Earth, which provides an 
angular resolution of 2.5°. The surface area covered on 
Earth is about 160,000 km^. The duty cycle of EUSO 
will be similar to that of surface fluorescent telescopes, 
of order of 10%. The plan is to have JEM-EUSO looking 
straight down and also in a tilted mode that will increase 
the viewing area by a factor up to 5 but decrease its res- 
olution. EUSO will also be equipped with devices that 
measure the transparency of the atmosphere and the ex- 
istence of clouds. Clouds are not always bad for such a 
detector because some shower signals could be reflected 
by them and detected better by the instrument. 

The motivation for EUSO is the study of cosmic rays 
of energy a bove 5x 10^^ eV as well as very high energy 
neutrinos ( (|Ebisuzaki 120081 )). 

Currently JEM-EUSO is the second experiment to be 
launched to the Japanese Experimental Module of ISS 
in 2015. The launch will be provided by the Japanese 
transfer vehicle HTV. 
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